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THE BOUNDARY-VALUE PROBLEM OF THE THEORY 
OF STELLAR ABSORPTION LINES 


BENGT STROMGREN 


ABSTRACT 


The problem of the theoretical calculation of intensities within absorption lines is con- 
sidered. The case when the ratio of continuous absorption to continuous-plus-line ab- 
sorption is variable through the atmosphere in such a way that the relative changes are 
not very large but are otherwise arbitrary is considered. For this case convenient formu- 
lae for the calculation of absorption-line intensities in the integrated spectrum are de- 
rived. With the aid of the relations derived, the relative influence of the various layers 
of the atmos phere on the intensities within absorption lines is discussed. The relations 
are specialized for comparison with equations previously derived by Eddington and 
Unséld. The accuracy and range of applicability of the procedure developed is discussed. 


I 


The theory of stellar absorption lines gives the relation between 
the ratio r of the intensity at any point within an absorption line to 
the intensity of the adjacent continuous spectrum, on one hand, 
and the ratio 7 of the line-absorption coefficient to the coefficient of 
continuous absorption, on the other. As is well known, the simple 
case that the ratio 7 is constant throughout the atmosphere can be 
treated analytically in a simple way. This case has been discussed 
by Eddington and Milne. Eddington’ has shown that the case, when 
1/(1 + ») is a linear function of the optical depth 7 in the continu- 
ous spectrum, can also be treated analytically, the ratio r being 
expressible with the aid of Bessel functions. Woolley? has used the 
method of computing 7, just referred to, in a number of cases of 
practical importance, and has also used the method in the case when 

'M.N., 89, 620, 1929. 2 Thid., 92, 482, 1932; 93, 691, 1933. 
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1/(1 +7) can be put equal to a linear function of the optical depth 
for the part of the atmosphere above a certain boundary depth, 
and equal to a linear function differing from the first below that 
boundary depth. 

When 7 is given numerically as a function of optical depth, the 
ratio r can always be found with the aid of numerical integrations. 
This method has been used by Pannekoek.’ 

In many problems of quantitative analysis of stellar absorption- 
line spectra it is of importance to be able to carry out calculations 
of r values in a convenient, and yet not too inaccurate, way. The 
simplest method is to choose the value of at some representative 
point—say at r = 1/3—and compute r as if » were constant 
throughout the atmosphere and equal to that value. This method, 
however, may lead to results that are considerably in error. Unséld4 
has developed a very convenient method for calculating accurate 
r values in the extreme wings of stellar absorption lines (r > 0.85). 
This method is of great value in the discussion of broad absorption 
lines for which contours can be measured. In the case of fainter 
lines, where only equivalent widths can be determined experimental- 
ly, the method naturally fails. Those parts of the contour for which 
r < 0.85 contribute largely to the equivalent width, except in the 
case of extremely weak lines, which in general would not be measur- 
able. 

Minnaert’ has developed Unsdéld’s method for the extreme wings, 
with the purpose of making possible accurate predictions of the vari- 
ation of the intensities in the extreme wings over the solar disk. 

Ten Bruggencate® has treated the general boundary-value prob- 
lem of stellar absorption lines in a purely analytical way, subdividing 
the atmosphere into a large number of separate layers. Wellmann’ 
has developed analytical methods for calculating r values and has 
shown that they can be used in practical cases. 

In the present paper we shall consider the case that the ratio 
1/(1 + 7) of continuous absorption to continuous-plus-line absorp- 


3 [bid., 91, 139, 519, 1930-1931; Pub. Astr. Inst. Amsterdam, No. 4, 1935. 
4Zs. f. Ap., 4, 339, 1932. 6 Tbid., 4, 159, 1932. 
5 Ibid., 12, 313, 1936. 7 [bid., 12, 140, 1936. 
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tion is variable through the atmosphere, in such a way that the 
relative changes of 1/(1 + 7) throughout the relevant parts of the 
atmosphere are not very large but are otherwise arbitrary. For this 
case convenient formulae for the calculation of r values are derived. 
It is shown that the changes of 7 can be rather large without serious 
error in the computed 7 value. Though cases may occur in which 
the variation of 7 is such that the method fails, a great many cases 
of practical importance are covered with sufficient accuracy. 

A number of examples of the use of the method here developed 
for calculating 7 values may be found in a forthcoming paper by 
the author on the hydrogen content of stellar atmospheres. 

The procedure here developed can be applied to the theory of 
the continuous spectrum for the general case when the absorption 
coefficient is an unrestricted function of depth and wave-length. An 
investigation of this problem is, however, outside the scope of the 
present paper. 

II 


The theory of stellar absorption lines*® leads to the following dif- 
ferential equation for the intensity J, within an absorption line: 


cos 6 = (kv + — k,B, — (1 el,OB, (1) 


Here /, is the intensity at a certain geometrical depth x below 
the surface of the atmosphere, in a direction making the angle 6 with 
the outward normal to the atmosphere (considered as plane-parallel). 
Further, p is the density, k, the coefficient of continuous absorption 
at the frequency », /, the coefficient of line absorption at frequency », 
B, the Planck intensity corresponding to the local temperature and 
J, the mean intensity (i.e., 7, averaged over all directions), both at 
frequency v. The quantity ¢ is a measure of the influence of electron 
capture on the upper stationary state of the line considered and of 
collision transition to that state upon the intensities within the 
line, while Q is a measure of the deviation from a Planck spectrum 

8 Cf. A. S. Eddington, The Internal Constitution of the Stars, Cambridge, 1926; E. A. 
Milne, Handbuch d. Astrophysik, 3, 1, Berlin, 1929; S. Rosseland, Astrophysik, 


Berlin, 1931; B. Strémgren, Handbuch d. Astrophysik, 7, Berlin, 1936; S. Rosseland, 
Theoretical Astrophysics, Oxford, 1936. 
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of the ionizing continuous spectrum at the point considered. The 
exact definition of the quantities « and Q can be found in a paper 
by the author;? in the present discussion we shall only make use 
of the fact that ¢ is a small quantity, usually less than 0.01, while Q 
is usually somewhat larger than 1 but of the order of magnitude r. 

The equation of transfer (1) is the most general equation for the 
case of lines originating from the ground state, when line fluorescence 
is neglected. 

When the system of differential equations (1) is treated according 
to Eddington’s method of approximation, the following equations 
are obtained. Averaging (1) over all directions, one finds 


dH, 
= (ky + Jy = (ky + , (2) 


H, being defined as 
dw 
H, = fl. (3) 
or 


H, = —F,, (4) 


where F, is the net flux. Averaging (1) after multiplication by cos 6, 
we get 


dK, 
(ky + (5) 
K, being defined as 
K, = cos (6) 


Introducing Eddington’s approximation K, = 3J, in (5), we obtain 


dJ, 
pdx 


= 3( + (7) 


9Zs. f. Ap., 10, 237, 1935. 
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THEORY OF STELLAR ABSORPTION LINES 5 


Now J, can be determined from (2) and (7), when use is made of 
the boundary conditions of the problem. It is convenient at this 
stage to introduce the optical depth in the line at the frequency v: 


or 
i= { (ky + pdx (9) 
The following simplifying notation is used: 
l, 
(10) 
= 
= 1+ . (11) 


Then equations (2) and (7) take the form 


dH, _ _itedQ 
dt, MJ, 1+ em B,) 
and 

dJy _ 

3H... (13) 


Eliminating H, between (12) and (13), we find 


The factor (1 + e7,Q)/(1 + en,), because of the smallness of e, 
does not differ appreciably from unity, except in the case of very 
large n, values, i.e., at frequencies near the center of strong absorp- 
tion lines. The results of the present investigation are intended 
mainly for use in calculations of equivalent widths of absorption 
lines. In these calculations the factor in question can be put equal 
to unity without appreciable loss of accuracy, the corresponding 
small changes of the intensity near the center of the lines leading to 


dty = (ky + 1») pdx (8) 
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very small changes in the equivalent width only. For the sake of con- 
venience we shall, therefore, in the following analysis omit the factor 
(1 + en.Q)/(1 + en,). If for other purposes it should be desirable 
to take account of this factor, one can do so simply by interpret- 
ing B, in the relations derived below as the Planck intensity cor- 
rected by multiplication with (1 + €7,Q)/(1 + e7,). We therefore 
write (14) as 
= 3.(J, — B,). (15) 
The proper particular solution of (15) is fixed by the following 
boundary conditions: (1) As ¢, — «, J, — B, must not increase ex- 
ponentially (in important simplified cases J, — B, then tends to 
zero). (2) For ¢, = o, J, = 2H, (Eddington’s approximate expres- 
sion of the fact that on the surface there is no ingoing radiation): 


(Jv) = 2(H»)1,-0- (16) 


Making use of (13), the second boundary condition can be written 
in the following form: 
= (17) 
When J, has been found as a function of ¢,, 7, can be determined 
as a function of ¢, and @ from (1). When it is sufficient to find the 
intensity in the integrated spectrum radiated by the star, it is not 
necessary to go through the latter analysis, as was pointed out by 
Eddington. The net flux F, for t, = 0, which is proportional to the 
intensity in the integrated spectrum, is, according to (4) and (13) 


or (16), given by 


and can thus be found immediately, when J, has been determined 
as a function of ¢,. 


| 
| 
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Hl 
The mathematical boundary-value problem of the theory of 
stellar absorption lines can thus be stated as follows: One has to 
find that solution of the differential equation 
3.(J> — B,) (19) 
which satisfies the boundary condition that J, — B, does not in- 
crease exponentially as ¢, goes to infinity, and the further boundary 
condition that 


The intensity within the spectral line in the integrated spectrum of 
the star is given by (J,),,-.. In (19) A, and B, should be considered 
as known functions of ¢,. 

For the sake of convenience we shall drop the index » of all 


quantities in the analysis following. Let 


y=J—-B. (21) 
Then, from (19), 
ay _ @B 
(22) 


The boundary conditions are now that y must not increase exponen- 
tially as ¢ goes to infinity, and, according to (20), 


Introduce 
(24) 
(25) 
then 
— = a — (26) 
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Finally, (18) becomes 


3 


Again \ and d’B/d? are to be considered as given functions of ¢. 
The simplified case first studied by Eddington and Milne corre- 
sponds to Xd equal to a constant, independent of ¢, and d?B/dt? = 0, 
i.e., B = a + pt. In that case it follows immediately that 


— 
3 


Then, from (27), 


or 


(30) 


The introduction of the values of a and p following from a dis- 
cussion of the Planck intensity as a function of the optical depth ¢ 
then completes the calculation. 

We shall consider in the following the case when the variation of 
\ is small compared with A. The boundary-value problem can then 
be solved with the aid of standard perturbation methods of the 
theory of differential equations. The first part of the analysis, up 
to equation (51), is, however, perfectly general. 

Let 


A=rA+6, (31) 


where X, is some properly chosen constant value, not necessarily the 
value of \ for t = o. Then (22) becomes 


_ 


+ = + 20. (27) 
| 
| 
| 
| 


THEORY OF STELLAR ABSORPTION LINES 9 
We can write the general solution of (32) in the form 
y = CeVit + (33) 
where C and D are functions of ¢. From (33) it follows that 


dy 


VB CoV! — V (34) 


provided we fix C and D by the usual condition 


Then 
v3 
di? — 3rot 3rA.De 3ro t 
(36) 


and therefore, according to (32), one must have 

V 3X0 V 3X. 


Solving (35) and (37) for the derivatives of C and D, one finds 


dC I = I d?B 
dt 2V 2V dt? (38) 
and 
dD I 


| 
| These two differential equations can be solved by quadratures: : 
| I I 
and | 
I I ~, @B 
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From (33), (40), and (41) one finds that the general solution of 
the differential equation (32) can be written in the form 


2V 3X. t dv 

(42) 
t 2 


The two constants of the general solution are c, and c,. The next 
step is to determine c, and c, in such a way that the boundary condi- 
tions are satisfied. As ¢ tends to infinity, the two terms involving c, 
will be of the order of magnitude of y. Hence, the sum of the two 
terms involving c, must not increase exponentially with ¢. Suppose 
that c, were finite. Then the sum of the two integrals involving c, 
would tend toward a finite limit as ¢ goes to infinity, and the sum 
of the two terms involving c, would increase exponentially as ¢ goes 
to infinity, in contradiction with the result derived from the bound- 
ary condition. Hence c, has to be replaced by ~. Then the first 
boundary condition is satisfied. 

Putting c, equal to zero, allowing for this by adding a constant 
times the factor of the integrals involving c,, one can now write the 
solution as 


2V 3X. t dt? 


> (43) 
3No 


t 2 
2VvV 3 ° 


where y is a constant of integration. 


| 
| 
i 
| | | 
| 
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We have now to determine y so that (43) satisfies the second 
boundary condition (26). From (43) one derives 


2 (dy I weer 
dt 


(44) 


3V 3X0 dt? 3V 3X0 
I 2 
and 
5 + dt 


| (45) 


Inserting (43) and (44) in the boundary condition (26), one gets 
the following equation for the determination of y: 


(46) 
— (1+ 3V3y)y= a — 3p. 


Solving this equation for y, one finds 


+ 2V 1 + 


2V 1 + Jo dt? 


a— 3p e~ V3! . dt 
° (47) 


Inserting (47) in (43), one arrives at the particular solution satis- 
fying the two boundary conditions. In order to find the quantity 
(F/7),-., Which is proportional to the intensity in the integrated 
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spectrum, one has to insert (47) in (44) or (45), and then (44) or (45) 
in (27). Carrying out the last-mentioned operations, one finds 


a — 3p Vn 


2V 3d I + 2V ° dv 


3r0t 
-6-dt 
val, 


dt 


F fav 3d. + 4p 4 { 
= -y-+d-dt 
( 1+ 20/2), 1+ 2/3), | 


4 
+ 2/2), ° dt 


Up to this point the analysis is perfectly general. It is seen that 
the boundary-value problem is solved when a solution y is known, 
accurate or approximate, which can be substituted into the first 
integral on the right-hand side of (49). 

In the simple case discussed on page 8, that \ is constant and 
d*B/dt? is zero, we find, on putting \ = A., that 6 is o. In that case 
(49) reduces to (30), as it should. 

Suppose now that 6 and d?B/d?? are small quantities of the first 
order. Then (28) with A = X, is an approximate solution y, the 
error of which is of the first order. Introducing this y in (49), one 
consequently finds the intensity in the integrated spectrum with an 
error that is only of the second order. 

In the same way, substituting the approximate solution (28) in 
(42), one gets the general solution y with an error of the second order. 
This is, in fact, the standard procedure for finding the perturbed 
general solution of (32). Similarly, the particular solution satisfying 
the boundary conditions could be found from (43) and (47), substi- 
tuting (28) in the integrals of these equations. 

If a better approximation to y than (28) is known, it can, of course, 


(48) 


1/F 
wad, dt? 
or 
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be used in (49) to give a still better approximation to the intensity. 
For instance, the Bessel function solution mentioned on page 1 
might be used for this purpose. 

Finally, it should be mentioned that if \ is constant, so that 
5 = o, the correct value of the intensity can be found from (49), ir- 
respective of the properties of the function B(#). This corresponds 

48) to the fact pointed out by Milne’ and Rosseland® that it is rela- 

tively simple to avoid assumptions concerning B in solving the 
boundary-value problem. 

The dependence of the intensity in the integrated spectrum upon 
B is perhaps more clearly brought out when the second integral in 
(49) is transformed by two partial integrations. One has identi- 
cally (cf. [24] and [25]): 


{ e~V3het — a — p+ anf e~V3et. Bs dt. (50) 


Inserting this identity in (49), we find 


-3-a 


This relation holds quite generally. In the special case that X is 
j constant, we have \ = \, and 6 = 0, so that the first integral 
| vanishes. The second integral is always proportional to a cer- 
tain mean value of B. The larger the ratio 7 of line absorption 
to continuous absorption, the smaller is \, and the larger is the 
exponent for the same optical depth 7 in the continuous spectrum 
[Vv tw (Vv 3. and the higher are the layers that contribute 
appreciably when forming the proper mean value of B. 
Substituting now the approximate solution (28) with \ = \, in 


% Phil. Trans. Royal Society, London, Ser. A, 223, 201, 1923. 
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(49), one gets the following expression for the intensity in the inte- 
grated spectrum, the error of which is of the second order: 


(52) 
+ — -dt. 
1+ 3V3dJo dt 

Similarly, (51) gives the equivalent expression 
F 

dt 

(53) 


In the form given, the equations are suited for numerical deter- 
minations of intensities within stellar absorption lines. For given 
values of \ and B as functions of ¢, one determines a and p, chooses 
a suitable \,, then finds 6 as a function of ¢ and carries out the 
quadratures numerically. For purposes of interpretation of the 
equations we shall, however, carry out a few further transformations. 


IV 


Let us consider for a moment the case that Xd is constant and 
d’?B/dt? is equal to zero. One might then choose a , differing from 
d by a quantity of the first order. The corresponding 6 would be 
constant and of the first order. Carrying out the quadrature for 6 
constant, we arrive at an expression for (F/7),-o, which must agree 
with (30), neglecting quantities of the second order in 6. These con- | 
siderations lead to the following transformations of (52). Introduce : 
the quantities 


q = (54) 
and 
| go = 3V3%, (55) 
and define Ag by 


q= q+ dq. 


| 
| 


| 
| 
| 
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Then, neglecting terms of the second order, 


26 
Ag =—. 
(57) 


The part of (52) which is proportional to a can now be re-written 
as follows: 


t+ 3V3r (1 + 


2a (58) 
=> —3qot . . 
Now, 
d/2qa\_ 2a 


Therefore one has, neglecting quantities of the second order, 


(1 + 3V 3%)? Jo 1 + + Aq)’ 


' where Ag has been defined by 
Defining @ by (cf. [56]) 
(62) 


we can consequently write the part of equation (52) that is propor- 
tional to a in the form 


(63) 
From (61) and (62) it follows that 


) 
| 
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Quite similarly one derives that the part of (52) proportional to p 
can be written as follows: 


* 


Inserting the expressions thus derived in (52), we get 


F\ _ 209 + 4p 


According to (66), the formal difference between the case of con- 
stant and that of varying ) is that in the latter case a certain mean 
value of V) has to be substituted for V.. It should be noted that 
it is only by introducing the mean value of a quantity proportional 
to the square root of d (cf. [54] and [64]) that one arrives at an equa- 
tion having this property. 

Assuming now in (66) that d?B/d?? is a small quantity of the first 
order, one can substitute g for go = 24/3), in the second term on 
the right-hand side of (66), the error being of the second order. Per- 
forming, further, two partial integrations (cf. [50]), we find 


F\ _ 209+ 4p 203+ 
or™ 


Introducing again X instead of g (cf. [54] and [55]), and inserting 
for the sake of clarity the indices v (cf. p. 7), the equations for 
calculating the intensity in the integrated spectrum are (cf. [64] 


and [68]) 


™ It is interesting to compare eq. (68) with an equation given by Milne (Joc. cit.,! eq. 
60]), which corresponds to the special case of \ constant throughout the atmosphere. 
The structure of the equations is exactly the same. 


| 

| 

| 

Wi 

| 
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and 


(71) 
y+ 


In the very simplest case that both \, and B, are constant, this 
reduces to 


B,, (dvand B, are constant), (72) 
as it should (cf. [30], with a = B and p = o). 

The procedure for calculating the intensity within the absorption 
line according to (69), (70), and (71) can be described as follows: 
For x, and B, varying with the optical depth t», we strike a mean of the 
values with the weighting function then a mean of the 
B, values with the weighting function e~V3V™'», The intensity is cal- 
culated as if X, and B, were constant throughout the atmosphere and 
equal to the mean values mentioned. 


V 


In many practical cases sufficient accuracy is obtained by putting 
B, equal to a linear function of the optical depth 7 in the continuous 
spectrum 


B, = dy + br (73) 


(rad 


k being the Rosseland mean of the continuous absorption coefficient 
over the spectrum. It should be noted that when \, is not constant, 
this may be an approximation of sufficient accuracy, while the ap- 
proximation of B, equal to a linear function of the optical depth ¢, 
within the spectral line (cf. [9]) is not. 

We shall consider in this section the form the equations developed 
for calculating the intensity within the absorption lines take when 
the assumption (73) is made. 


where 


| = 
| 

q 
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From (9) and (74) it follows that 


ad k+l, 
(75) 
Defining 1, as the ratio of k, and k, 
ky = nok (76) 
we have (cf. [10]) 
dt, 
n(1t + m). (77) 


In the analysis of the continuous spectrum it will usually have 
been assumed that , is constant throughout the atmosphere. We 
shall, for simplicity, make the same assumption here, though it pre- 
sents no difficulty to treat the case when n, is variable. 

From (73) and (77) it follows that 


dB, _ 


(78) 
Therefore, 
B, = a + Ny I+” dt, (79) 
Call 
(80) 
so that 
B, = dy + pt, (81) 


In the case that € = o, yw, thus defined coincides with \,. The 
difference between the two quantities generally will only be appre- 
ciable for large values of 7,, i.e., near the center of strong absorption 
lines. 

We now introduce the expression (81) for B, in equation (70), 
defining B,, 


B, = d + Vy f dt, . (82) 


a 
i 
I, 


ALES 
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The second term on the right-hand side of (82) is transformed by a 
partial integration: 


— 
B=a+t+ + dty. (83) 


Jo 


Defining a mean value of yu, by 


we have 
(8s) 


Introducing this expression for B, in (71), one finally gets 


F avr, +4 by 
/ t,=0 I be 2/3 Vx, 


In the case considered in this section, that B, = a, + 6,7, we 
therefore find that the intensity within the absorption line is calculated 
by striking proper averages (cf. [69] and [84]) over the quantities 
= (1 + ey,)/(1 + 9.) and p, = 1/(1 + and then introducing 
these into the simple formula valid when X, and p, are constant through- 
out the atmos phere. 

In calculating equivalent widths of stellar absorption lines € can 
usually be put equal to zero. The corresponding small change of 
the contour near the centerof the line means only a small change in 
the equivalent width. In that case (86) reduces to 


4b, — 


F, 
= B, = by = 8 
Vd, ( ay + br € ( 7) 


with VX, and }, given by (69) and (84). 


| 
t 
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It means neg'ecting quantities of the second order only, if in (69) 
and (84) we put 


I 
V(dv)o ly = (88) 
and 
Therefore 
Vin = (Vie a( 5) (00) 
and 
V3 
— _ ) 
{ a( (01) 


In (87) one can substitute the usual expressions® for a, and b, ob- 
tained by expanding B, in a Taylor series from the surface, assuming 
the temperature to be given by 


= T4(1 + 47), (92) 


where 7, is the surface temperature. 
The contrast 
F, 


r= Freon) (93) 
between spectral line and continuous spectrum is obtained by put- 
ting VX, and X, in (87) equal to 1 in order to get FS". One thus 


finally gets 


$V34+ I+ 3V3 
where 
hy 
Xo = a (95) 


q 


| 
rv 
i 
| 
| 
| 
| 
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For many purposes, equations (go), (91), (94), and (95) will be 
convenient and sufficiently accurate for calculating r,. 

It is interesting to note the difference in the weighting functions 
for VX, and A, (cf. [90] and [91]). Going from the surface inward, 
the weighting function for VX, falls off much more rapidly than 
that for \,. This means that for ‘wo reasons the intensity near the 
center of an absorption line depends upon the properties of layers 
higher than those that determine the intensity in the wing of the 
line. First, as we go from the wing toward the center, V(d,)o de- 
creases and the weighting factors for both VX, and X, fall off more 
and more rapidly with increasing 7. Second, as one approaches the 
center, the \, term in (94) becomes smaller and relatively less im- 
portant compared with the V\, term. Therefore, the effective 
weighting function tends toward that for VA, as we approach the 
center, and the latter weighting function is the one that falls off 
more rapidly. 

The difference between the weighting functions just referred to 
is perhaps best interpreted by considering the expression (71) for 
the intensity within the absorption line. Quite general y, line ab- 
sorption affects the intensity within the radiated spectrum in two 
ways. First, it cuts down the intensity as would increased continu- 
ous absorption, the effect depending upon the existence of a tempera- 
ture gradient in the stellar atmosphere. Second, it causes a devia- 
tion from Kirchhoff’s law, in the sense that emission is reduced, so 
that there is a further reduction in the radiated intensity. The first 
effect can be studied formally by putting « = 1 (cf. the equation 
of transfer [1]) and hence \, = 1 (cf. [11]). Considering (70) and 
(71), one now sees that the first effect is associated with the weight- 
ing factor that falls off relatively slowly, whereas the second effect 
is connected with the weighting factor that falls off rapidly. 

The first effect, that analogous to the effect of large continuous 
absorption, cannot reduce the intensity below the Planck intensity 
corresponding to the surface temperature (cf. [70] and [71]). There- 
fore, when near the center of a strong line the line absorption is 
already great, a change in the line absorption means relatively little 
as far as the first effect is concerned. This corresponds to the de- 


i 
| 
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creasing importance of the \, term in (97) as one approaches the 
center of the absorption line. 

Finally, it should be mentioned that the first effect is naturally 
the more important the greater the change of the Planck intensity B, 
with optical depth 7, and therefore is more important the larger the 
ratio hv/kT, (cf. the increase of the limb darkening with decreasing 
wave-length). Corresponding to this, x)» occurs as a factor of the 
\, term in (94). The factor 1/n, in the same term can be interpreted 
in a similar way. 


VI 


The case of the extreme wings of absorption lines studied by 
Unsédld is a special case of that considered in the present paper. For 
in this case 7, is a small quantity of the first order, and hence X, dif- 
fers from a constant, namely unity, only by a small quantity of the 
first order throughout the atmosphere. 

In this section we shall specialize the equations obtained in the 
previous sections to the case that m, is small. The comparison with 
Unséld’s result is of interest, because Unséld did not make use of 
the Eddington approximations in treating the infinite system of dif- 
ferential equations (1) but considered the corresponding Schwarz- 
schild integral equation. 

We shall limit ourselves to the most important case when € = o 
and B, = a, + b,7. We can then determine 7, from (go), (91), and 
(94), putting A, = 1 — », and (A,). = 1. Thus, one finds 


I Xo ) 
(r+ 3V3)(3V3 
4n 
> (96) 
-m:dr. 
2 
3 | 


| 
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Writing Unséld’s formula in the form 


= Bia 3 Xo Bis) dr 
I Xo Jo 8 n 
I + v0 Jo 4 
4 Nn 


it is easy to see the correspondence of the various terms. The ex- 
ponential integrals correspond to e~¥3* and the square or the prod- 
uct of the exponential integrals corresponds to e~?V3*. A similar cor- 
respondence was noted by Minnaert in the investigation mentioned 
on page 2. In this connection it is interesting to note that in the 
theory of the continuous spectrum of a gray body, Eddington’s ap- 
proximation leads to an integral expression for the flux in terms of 
the Planck intensity, which differs from the Schwarzschild-Milne ex- 
pression® only by the substitution of e~¥%* for Ei,(r). 

Numerical values of the weighting function G for y, according to 
(96) and (97) are given in Table 1 for x,/m equal to 2, 4, and 6, re- 
spectively. It is seen that the differences between the functions are 
small. 


TABLE 1* 

n n n 

G G G G G G 

= Acc. (96) | Acc. (97) " Acc. (96) | Acc. (97) * Acc. (96) | Acc. (97) 

° 0.80 0.92 ° 0.80 0.88 ° 0.85 
o.1 o.1 68 .66 o.1 .70 .69 
0.2 .47 0.2 0.2 .60 .58 
0.5 .26 0.5 34 33 6.5 .38 .38 
1.0 .10 .12 1.0 .14 -16 1.0 
.040 .058 .060 .085 .072 .10 
2.0 .O17 .039 2.0 .025 2:0 .030 .054 
3.0 | 0.0029 | 0.009 3.0 | 0.0045 | 0.013 3.0 | 0.0054 | 0.016 


* Weight function for the extreme wing: moe f G-ny,-dr. 
° 


24 BENGT STROMGREN 


VII 
As was mentioned on page 1, Eddington has treated the case 


when 1/(1 + 7,) is a linear function of the optical depth r, 
a + Br, (98) 


with 
B, = dp + br (99) 


We shall limit ourselves again to the case that « = o. Also we shall 
consider the case only that x,/n = 4, so that b, = 3a,. Then the 
expression for 7, is' (cf. also Woolley?) 


= Ko(V120/6") (100) 
2 Ki(v1 2a/?) 
1 — 3V3a- 
K,(V12a/B?) 
Here” K, is a Bessel function of order zero and K} = —K, is the 


derivative of this function. Tables of these functions are given in 
Watson’s Bessel Functions (pp. 698 and 737) and also in Jahnke- 
Emde’s Tables of Functions (2d rev. ed., p. 286). (In the notation 
of Jahnke-Emde, they are proportional to iH (iV 1120/8?) and 
—H (iV 12/6"), the functions H® and H being Hankel’s func- 
tions of the third kind, of orders 0 and 1.) 

For B = 0, equation (100) reduces to (94) with x,/n = 4 and 
with \, constant throughout the atmosphere. The case X, = 
1/(1 + »,) given by (98), with 6 small, can be treated with the aid 
of both (94) and (100). We shall now consider this case and shall 
show that the results agree. This means that we specialize (94) for 
the case when the variation of \, assumed small but otherwise arbi- 
trary, is linear, and specialize (100) for the case when the variation 
of X, assumed linear but otherwise unrestricted, is small. 


12 Eddington considers in detail only the case when 8 is positive. In this case eq. 
(100) holds true. Thus the following derivation of (103) is true for 8 positive. For 8 
small it is very easy to see that in the case of negative 6 the Bessel functions J, and J; 
have to be substituted for Ky and K; in (100). Further, in this case J;/J, is, apart from 
terms of the second order, equal to the reciprocal of K:/Ko. It follows that (103) holds 
true for negative values of 8 also. 
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When 8 is small, the argument of the Bessel functions is large. 
Then, as was pointed out by Eddington, we have 


AV 12a/B? i2a 
K,(V120/B") I ( > :) . (101) 


8V12a/? 
In the limiting case considered, we may write this as 


_ 


4 Noa’ 


Introducing this expression in (100) and neglecting quantities of 
the second order and higher, we find 


= (103) 


Ny = 


On the other hand, on introducing (98), which, neglecting quan- 
tities of the second order and higher, is equivalent to 


a+ Bt, (104) 
in (go) and (91), we find 
Vy = Vv . 
and 
X= a+ (106) 


r= (107) 


which is identical with (103). 


24/ 
I+ 3 a 4 V 
Inserting this expression in (94), we obtain, for x,./n = 4, 
| V30(1+ 4 0+ 
| 
| | 
| 
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It is interesting to note that (105) can be written in the form 


Vy, a + B (108) 


which is equal to the square root of the value of d, for 7 = V a/2V ‘. 
Similarly \,, according to (106), is equal to the value of \, for r = 
Va/V3. 

It follows that in the case of a small linear variation of X, with 7, 
the proper mean values to use in (94) are: for V4 ,, the square root of 
the value of », at r= V (ry)o/2V 3; and for X,, the value of X, at 
V 3. 


TABLE 2* 
Ve} = r 
= 4= r 

0.56 0.28 0.88 0.96 0.94 0.04 
52 26 65 .86 82 82 
0.36 0.18 0.25 0.56 0.54 0.55 


‘ — values have been slightly revised, using J:/Jo instead of K:/Ko because 8 is negative 

The discussion on page 21 shows that in this case the representa- 
tive point moves from toward 3 as we pass 
from the wing of the line toward the center, there being thus two 
causes for the shift toward higher layers, namely, the effect just 
mentioned, and the effect of the decrease of V(X,)o, which quan- 
tity is roughly equal to 7,. 

Table 2 illustrates the use of the rule just stated. The examples 
are taken from an investigation by Woolley’ and refer to the case 
of the Ca lines in the solar spectrum. 

In all these examples there is agreement practically within two- 
figure accuracy. 

It is interesting to determine the accuracy of 7, values computed 
from (go), (91), and (94) in cases of rather large linear variations 
of X,, making use of (100). Table 3 gives two typical examples. In 
spite of the very considerable variation of 7, with depth, the results 


= 
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are only slightly in error. This is connected with the fact that X, 
has the limits o and 1. 


TABLE 3 
r ny at ny at 
a v v v 
\ B = vA, Acc. (94) T=0 T=0.5 
I 0.35 | 0.18 99 4.5 


The relative change of A, from the surface to the representative 
point 7 = V (d,)o/ 2V 3 is a measure of the variation of X, in the 
relevant part of the atmosphere: 


Relative change = * - = 
2V3 


(109) 

Comparing with (100), we see that this quantity is the reciprocal 
of the argument of the Bessel function in the expression for 7,. Thus, 
when the relative change is small, the argument is large and the 
simplified expression (101) is valid. 

Quite generally one gets a good estimate of the accuracy of r, 
values obtained with the aid of (90), (91), and (94) by calculating 7, 
first with (A,)o equal to the value of A, on the surface, and then 
with (A,)o equal to (Vv r,) . As was stated in the introduction, the 
accuracy proves to be sufficient in a number of cases of practical 
importance. 


YERKES OBSERVATORY 
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REGIONAL STUDY OF INTERSTELLAR SODIUM* 
PAUL W. MERRILL 


ABSTRACT 


Introduction.—The main circumstances which underlie an investigation of inter- 


stellar lines are briefly discussed. 

Line intensity and distance.—Data for double stars (Table 1) bear out the expecta- 
tion that interstellar lines should behave alike in the spectra of the components of a 
pair. The relationship between distance and intensity of the D lines was examined in 
23 small groups of stars scattered along the Milky Way from longitude 310° through 
0° to 190°. In several groups the more distant stars exhibit little increase in line in- 
tensity. Possible causes of this effect are discussed. Normal places from the 23 areas, 
when brought together (Fig. 2), show that regional differences, while probably present, 
are not extremely great. The mean relationship between intensity, y = 0.5(D2 + Dr) 
and x = distance corrected for space absorption, is represented to 2000 parsecs by 


the equation 
y = 0.001440 2° 8 | 


The distance-intensity relationship in seven large regions (Figs. 3 and 4), chosen with 
regard to the null points of galactic rotation, indicates that the standard galactic rota- 
tion plays a relatively minor role in producing the observed intensities of detached D 
lines. Little evidence of an increase of density of sodium vapor toward the galactic 
center appears in the region observed. Stars in areas with strong reddening do not be- 


have abnormally. 

Line intensity and displacement.—Longitudes 311° through 0° to 212° were divided 
into 14 regions, for each of which a correlation diagram between residual radial velocity 
and intensity is exhibited in Figures 5 and 6. In spite of considerable scatter, the in- 
fluence of galactic rotation is plainly marked in the observed motions. 


INTRODUCTION 

Numerous investigations have demonstrated a statistical increase 
in the intensities of detached lines with the distances of the stars in 
whose spectra the lines appear, and this relationship has been used 
to find the distances of early-type stars and of novae. While the 
method appears to be fairly successful, its reliability for individual 
objects remains to be determined. For this reason, as well as to in- 
crease our general knowledge of the properties of interstellar gas, a 
detailed investigation of the relationship between distance and line 
intensity is desirable. 

The intensity of an interstellar line, say D2 of sodium, in the 
spectrum of any star depends not only on the total number of sodium 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 


ington, No. 569. 
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atoms (in the ground state) along the line of sight but also on their 
motions. Data on displacements should therefore be included in a 
general study of interstellar lines. Direction, however, as well as dis- 
tance, must be considered because there is no guaranty that densities 
and velocities have a spherically symmetrical distribution about the 
sun; and effects varying in some systematic way with direction, as 
well as others due to spatial irregularities in density or motion, may 
be present. 

Studies of galactic rotation suggest that, for a given distance, in- 
tensities of interstellar lines may be some function of 


sin 2(/ — 331°), 


where / is the galactic longitude of the light path, 331° that of the 
center of rotation.’ In addition, since distant stars in high galactic 
latitudes probably lie outside the gaseous stratum along the plane 
of the galaxy, the intensities in such stars should be related to the 
cosecant of the latitude. This latitude effect appears to be impor- 
tant, however, in relatively few of the stars on the Mount Wilson 
program, which means simply that most of the brighter B-type stars 
lie within the boundaries of the gaseous stratum. The emphasis in 
the present investigation is therefore placed on the behavior of in- 
terstellar lines with respect to longitude. Distant stars in high lati- 
tudes should be discussed separately. 

Irregularities may be caused by clouds of interstellar gas having 
a partly random distribution, perhaps analogous to that of the dark 
dust clouds that obstruct our view of the Milky Way; or by local 
stream motions comparable with those of relatively small groups of 
stars; or by differences in kinetic motions of atoms in various por- 
tions of space. 

The observational investigation of these problems is not alto- 
gether simple. The ideal set of data would include intensities and 
displacements of interstellar lines in the spectra of a large number of 
stars uniformly distributed throughout the galaxy. Practical con- 
siderations, however, greatly restrict the number of suitable stars. 
In the first place, a serious limitation arises from the small distances 


"A. S. Eddington, M.N., 95, 11, 1934. 
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of most of the brighter stars relative to the dimensions of the galaxy. 
Among objects brighter than apparent magnitude 7.0, only the 
scarce c stars have distances greater than 2000 parsecs. Moreover, 
stars whose own spectral lines prevent reliable measurements of the 
interstellar lines must be excluded. This requirement eliminates 
practically all stars later than B3 for observations of detached H 
and K, and those later than Ao for the D lines, although somewhat 
later types may be of service in special circumstances. In the spectra 
of the distant c stars, for example, useful measurements of detached 
sodium lines can be made in types as late as A2. 

A practical difficulty in the discussion of the data is introduced by 
the fact that the distances of individual early-type stars are not 
known with accuracy. Distances derived from the differences be- 
tween the apparent magnitudes and the assumed absolute magni- 
tudes serve well for groups of stars but are not wholly reliable for 
individual objects. 


LINE INTENSITY AND DISTANCE 


A convenient first test of the behavior of interstellar lines is af- 
forded by observations of double stars. The components of a given 
pair are at practically the same distance from the observer, while 
the two light-paths are separated by very small angles. Correspond- 
ing interstellar lines in the spectra of the components would there- 
fore be expected to have equal intensity, and all interstellar lines 
should have the same displacement. With possible minor exceptions, 
these expectations are borne out by the measurements recorded in 
Table 1, for which the differences between corresponding values are 
not in excess of errors of measurement. The results for calcium de- 
pend on measures published by the Dominion Astrophysical Ob- 
servatory and on additional data obtained by Mr. R. F. Sanford 
and the writer at Mount Wilson, while those from sodium are wholly 
from observations at Mount Wilson. In certain pairs the velocities 
from the D lines are algebraically larger than those from H and K. 
This tendency, previously noticed in a larger collection of data,? may 
be a regional effect. 

As already remarked, the main purpose of the present article is 


2 Mt. W. Contr., No. 564; Ap. J., 85, 73, 1937- 
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the comparison of intensities and displacements of detached sodium 
lines in the spectra of stars in various galactic longitudes. The in- 
tensities were determined by Mr. O. C. Wilson and Miss Cora G. 
Burwell from spectrograms taken by the writer with a plane-grating 
spectrograph at the Cassegrain focus of the 1oo-inch telescope. To 
make the results homogeneous, that part of the process where per- 


TABLE 1 
DETACHED LINES IN SPECTRA OF DOUBLE STARS 
INTENSITY VELOCITY 
(EA) (Ku /SEc) 
STAR HD ADS | Sep. | Mac. | Type 
D2 Di Ca Na 
28446 Br | 3274} 1073} 5.4] Bin | 0.36 | 0.34 |— 3.0/— 0.9 
BG 6.2 | Bo 36 |+ 4.4/+ 0.1 
36861 4179] 4.4| 3-7 | O08 42 .24 |+13.3/+19.9 
5.6 | Br .40 | .23 |4+15.3/-+22.1 
Oric..... 135 5.4 | OF 26 .16 |+20.7/+22.0 
37041 4188} 52.6) 5.2] Oo 32 |+13.1|+22.2 
6.5 | Br 26 12 |+11.2/-+15 
37742 4263 2.6} 2.0] Bon .17 |+14.0|+18.0 
BISCO 144217 9913} 13.7] 2.9 | Br .20 .12 |—13.6|—12.8 
147933 10049} 3.4] 5.2 | B3n .20 — 7.9 
+35°3930...| 190429 Br | 13312] 2.0] 7.2] |......]...... —10.9/— 6 
214168 16095} 22.4, 5.8] Bgne| .26] .18 |—11.3/—14.9 
6.6 | B2 0.24 | 0.16 |—12.3/—13.8 


sonal equation might enter was done entirely by Mr. Wilson. Most 
of the displacement measurements were made by Miss Burwell. 
The results for individual stars, together with a description of the 
methods employed, will be given in another Contribution. 

To ascertain whether the relationship between distance and in- 
tensity is subject to marked regional effects, 23 fairly small groups 
of stars were selected for study. These groups are scattered irregu- 
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larly along the Milky Way from longitude 310° through o° to 190°; 
most of them are rather near the galactic plane, only 5 having cen- 
ters with latitudes exceeding 6°. The distances were computed from 
the formula 


log D = 0.2(m— M) +1, 


where the apparent magnitude, m, was taken from the Henry Draper 
Catalogue, and the absolute magnitude, M, from Table 2, which 
represents a smoothed mean of the fairly consistent values adopted 
by several recent investigators. 


TABLE 2 
ASSUMED ABSOLUTE MAGNITUDES 
Lines Diffuse 

Type Others 
ess —3.4 —3.4 
—3.1 3-3 
—5.0 


For each of the 23 groups a correlation diagram was made between 
distance and the arithmetical mean of the measured intensities of 
D2 and Dr. Figure 1 shows sample diagrams for the following 
regions: h and x clusters in Perseus, longitude 102°; x Aurigae, 139°; 
Trifid Nebula, 337°. 

In several of the curves, the portion corresponding to the greater 
distances is nearly flat. In the x Aurigae region, for example, the in- 
tensities show but little increase beyond 800 parsecs. Possible causes 
of such an effect are: 

a) Obscuration.—If an absorbing medium is present, the calcu- 
lated distances of the more distant stars will be too great relative to 
the nearer ones, thus flattening the curve. For a uniform absorption 
of 0.35 visual magnitudes per kiloparsec, the displacement of the 
curve is not very marked for stars nearer than 1500 parsecs. Accord- 


| 
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i ing to Stebbins and Huffer,’ stars in the region of x Aurigae are sub- 
ject to considerable reddening, the average color excess of those 
plotted in Figure 1 being about 0.2 mag. (Stebbins’ C, scale). If the 
change of the visual magnitude by interstellar absorption is nu- 
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Fic, 1.—Relationship between distance and D-line intensity for individual stars in 
} three regions. 


merically equal to the accompanying color excess, as certain photo- 
metric investigations indicate, the shape of the observed curve con- 
necting distance and D-line intensity would not be greatly altered by 
| the correction for absorption. It seems improbable, therefore, that 
| the flattening is caused wholly by homogeneous absorption. 


3 Pub. Washburn Obs. Univ. Wisconsin, 15, 217, 1934. 
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b) The curve of growth of line intensity with number of atoms may, 
under certain circumstances, have a flat portion. Discussion of this 
point is reserved for another Contribution. 


TABLE 3 
INTENSITIES OF INTERSTELLAR SODIUM ‘ 
LINES IN VARIOUS REGIONS | 
No. of Distance Int. EA 
Stars (Parsecs) 0.5 (D2+D1) 
— 4° 4 600 0.23 
+ 1 II 1800 0.74 
+ 1 4 1100 0.55 
° 7 300 0.15 
1200 0.60 
—17 3 600 0.24 
1 II 700 0.51 
1700 0.77 | 
2600 
° 7 3600 0.84 ‘ 
— 3 18 1400 0.65 
2200 
3200 
+ 2 7 1200 0.64 
° 7 700 0.38 
1100 0.46 
+ 2 9 500 0.30 
1700 0.48 
—12 6 400 0.25 
—17 9 600 0.22 
+ 1 8 800 0.45 
+ 1 6 800 0.30 
fo) 6 600 0.35 
+19 8 200 0.17 
— 2 8 600 0.25 
1100 0.39 
+11 4 550 0.37 


c) Dispersion in the absolute magnitudes for a given type, or errors 
in the assumed relative magnitudes of various types, would cause 
stars actually at the same distance to be assigned a considerable 
range in distance, and the curve might thus be given a fictitious 
flattening. In other words, the regression-curve of intensity against 
distance may be less steep than the true relationship. 


Equivalent width 0.5 (D2+D1) 
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d) Irregularities in density of sodium vapor.—If a number of stars, 
situated in a relatively tenuous region of space, were seen through a 
sodium cloud of high density, the intensities would not progress with 
distance, and the plotted curve would be flat. This possibility should 
be given serious consideration. 

Since in most of the 23 regions the observations do not define 
clearly the way in which line intensity increases with distance, a 
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fo) 800 1600 2400 3200 parsecs 


Fic. 2.—Relationship between distance and D-line intensity for mean values derived 
from 23 small regions. The galactic longitude corresponding to each point is given. 
The continuous line is the estimated mean curve as observed; the dashed line, that cor- 
rected for space absorption of 0.35 mag. per kiloparsec. 


single normal place was formed for each whose co-ordinates are 
practically the means of those of the individual stars. Where the 
progression of intensity with distance is clearly exhibited, however, 
two or three pairs of values were determined. All these normal places 
are listed in Table 3 and plotted in Figure 2. 

The correlation between distance and intensity is sufficient to 
show that regional differences, while probably present, are not ex- 
tremely great. A single mean curve, such as that represented by the 
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continuous line, used for determining the distances of stars near the 
galactic plane, should give results of the right order of magnitude, 
although for a few areas better values may perhaps be obtained 
from the original regional curves. 

The dashed line shows the position of the curve when corrected 
for a space absorption of 0.35 mag. per kiloparsec. According to this 
curve, the intensity increases nearly, but not quite, in proportion to 
the distance. Up to 2000 parsecs a good representation is given by 
the equation 

y = 0.0014402x° 


x = distance in parsecs , 


y 


intensity, 0.5 (D2 + D1), in 
equivalent angstroms. 


Since the 23 small areas yield reasonably consistent results for the 
distance-intensity relationship, all the stars observed were next 
grouped according to longitude in seven larger regions chosen with 
regard to the maxima and the null points of galactic rotation. Three 
regions, each extending over 20° of longitude, are centered on the null 
points observable from Mount Wilson, namely, those at 331°, 61°, 
and 151°; two regions 70° wide are centered on the maximum points, 
16° and 106°; while the terminal regions, 307°-320° and 162°-212°, 
complete the list. Plots of the correlation between distance and D- 
line intensity are shown in Figures 3 and 4. 

Because the D lines are intrinsically very narrow, high total in- 
tensities can be produced by a moderate number of atoms only if 
some external effect causes the absorption to be spread over a range 
of wave-length greater than the normal. Such changes in wave- 
length may be caused by the influence of near-by particles on the 
absorbing atom, but this effect is probably negligible at the low 
densities of interstellar space. Gross electric or magnetic fields can 
scarcely be postulated. The only other known cause of widening is 
the Doppler effect arising from the relative motions of atoms and ob- 
server. It thus appears that differential motions along the line of 
sight are of essential importance in producing the observed intensi- 
ties of interstellar lines. 

If the predominating motions are those corresponding to general 
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galactic rotation of a planetary type, the interstellar lines should, for 
a given distance, be less intense in stars near the nodes or null points 
at 331°, 61°, 151°, and 241° than at intermediate longitudes. For dis- 
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Fic. 3.—Relationship between distance and D-line intensity for individual stars 
grouped according to galactic longitude. An arc above a plotted point indicates that 
the star is 10° or more north of the galactic equator; below, 10° or more south. The mean 
curve (observed) from Figure 2 is placed on each diagram for reference. 

Region 321°—341°.—The point at the upper right-hand corner of the diagram corre- 
sponds to HD 160529, cAgea. The sodium lines appear double, the probable explana- 
tion being that the component toward the red is interstellar, while that toward the 
violet, which is less intense, is caused by stellar absorption displaced as in P Cygni. 
Since the wings of the components overlap, some uncertainty exists as to the intensity 
of the interstellar line. The plotted value probably represents an upper limit, the true 
intensity being perhaps 20 per cent less. 

Region 342°-0°-50°.—The point at x = 2090, y = 0.37, corresponds to HD 190066, 
cBo. If it is an ordinary Bo star instead of a c star, its calculated distance would be only 
g6o0 parsecs. The points at x = 1580, y = 0.37, corresponding to HD 193077, Oc, and at 
x = 1900, y = 0.40, HD 191566, B2, have low weight. HD 227836 yields a point beyond 
the limits of the diagram at x = 5250, y = 1.43. 


tances up to 600 parsecs, two or three small areas near null points 
yield slightly low intensities (Fig. 2), but the data have low weight 
and the plots in Figures 3 and 4 make the effect appear uncertain. 
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At greater distances, where the effects of rotation should manifest 
themselves more strongly, areas at 337° and 156° (toward the galac- 
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Fic. 4.—Relationship between distance and D-line intensity for individual stars 
grouped according to galactic longitude. An arc above a plotted point indicates that the 
star is 10° or more north of the galactic equator; below, 10° or more south. The mean 
curve (observed) from Figure 2 is placed on each diagram for reference. 

Region 72°-140°.—The point at x = 1320, y = 1.27, corresponds to HD 223385 
(Boss 6111), cA2ea, which may be at a greater distance than that assigned to it. 
HD 4717 yields a point beyond the limits of the diagram at x = 5750, y = 0.08. 

Region 162°-212°.—The point at x = 1900, y = 0.39, corresponds to HD 46106, 
Bo (Victoria). The calculated distance for type B8 (Harvard) is 830 parsecs. 


tic center and in the opposite direction) exhibit low intensities 
(Fig. 2), while an area at 55° (near the intermediate null point) does 
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not. Figures 3 and 4, based on slightly more extensive data, display 
the same behavior. 

These results appear to indicate that the standard galactic rota- 
tion plays a relatively minor role in producing the observed inten- 
sities of detached D lines. This conclusion does not mean that inter- 
stellar matter fails to take part in galactic rotation or that the ex- 
pected widening of the lines is nonexistent. The interpretation is 
rather that the effect is badly blurred by interstellar motions of 
other kinds. 

A curious irregularity in which galactic rotation can scarcely play 
any part is the low intensity of the lines in the region of longitude 91° 
compared with those in regions 84° and 102°. Another interesting 
feature is the very low intensity found in certain stars at high galac- 
tic latitudes, in the region of longitudes 162° to 212°. This indicates 
either marked local irregularities or a surprisingly thin galactic stra- 
tum of sodium vapor. 

There is little evidence of a general increase of density of sodium 
vapor toward the galactic center; but the distances to which our ob- 
servations extend are, of course, small compared with the galactic 
radius. Moreover, it is probable that the intense lines in the more 
distant stars of our program would not be very sensitive to differ- 
ences in interstellar density. 

Points corresponding to areas in which Stebbins and Huffer have 
noted strong reddening’ are not distinguished by large residuals in 
the distance-intensity relationship. Intensity does show a positive 
correlation with color, as found by E. G. Williams for the K line;4 
but this relationship, if arising statistically from the fact that both 
intensity and color increase with distance, does not necessarily indi- 
cate a physical dependence. 

A general discussion of the distance-intensity relationship is 
planned for another Contribution. 


NOTES ON TWO STARS IN NEBULAE 


1. HD 164492 in the Trifid Nebula, M 20, NGC 6514.—The mean 
intensity of the D lines in the spectrum of this star in the southern 
lobe of the Trifid nebula is 0.16 equivalent angstroms, a value de- 


4 Mt. W. Contr., No. 487; Ap. J., 79, 280, 1934. 
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pending on one spectrogram of fair quality. The distance read from 
the curve in Figure 1 is 380 parsecs, which may be taken also as a 
determination of the distance of the nebula. With this distance and 
the assumed absolute magnitude of — 3.4 for spectral type O7, the 
apparent visual magnitude should be 4.5. The value in the Henry 
Draper Catalogue is 6.9, indicating an obscuration of 2.4 magnitudes. 
This result is not considered highly accurate but is given as an il- 
lustration of an interesting method. 

2. HD 164794 in M 8, NGC 6523.—For this star (which is nof 
the nuclear star) the intensity 0.5 (D2 + D1) = 0.31 EA; distance 
800 parsecs; type O5; M, —3.4. From these data the calculated ap- 
parent magnitude is 6.1, which agrees within errors with the HD 
value, 5.9. M 8 lies within a degree and a half of the Trifid Nebula. 


LINE INTENSITY AND DISPLACEMENT 


Doppler effects that are symmetrically distributed about the nor- 
mal position of the line center will increase the total line intensity 
without producing any displacement of the line as a whole. This 
would be true of effects due to kinetic motions or to other types of 
random motion whose projections on the line of sight have similar 
distributions in both positive and negative directions. It would not 
be true of motions such as differential galactic rotation whose pro- 
jections are wholly or predominantly of one sign; these motions 
give rise to displacements having a systematic relationship to in- 
tensity. 

A direct comparison of observed intensities and displacements at 
various longitudes is therefore in order. For this purpose the longi- 
tudes from 311° to 212° were divided into 14 regions, for each of 
which a correlation diagram for residual radial velocity and total in- 
tensity is exhibited in Figures 5 and 6. The residual velocity was 
computed by correcting the observed velocity for solar motion of 
20 km/sec toward R.A. 18"o™, Dec. +30°. 

It is instructive to compare the diagrams in sequence of longitude, 
bearing in mind that effects of galactic rotation should vanish at the 
null points, 331°, 61°, and 151°, and should have positive maxima at 
16° and 196° and a negative maximum at 106°. Actually for regions 
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which include null points there is little progression of displacement 
with increasing intensity; for regions near the maxima, the increas- 
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Fic. 5.—Relationship between residual radial velocity and intensity of the D lines 
for individual stars grouped according to galactic longitude. An arc above a plotted 
point indicates that the star is 10° or more north of the galactic equator; below, 10° or 
more south. The dashed line on each diagram marks the locus at which displacement 
equals one-half the equivalent width. 

Region 34°-44°.—Many of the stars included in this area are faint, and the data have 
low weight. Nevertheless, the plot is considered representative. The very divergent 
point « = 0.14, y = +24, corresponds to HD 178849, B3, mag. 6.6. The data, derived 
from one good spectrogram, need confirmation. 


ing displacements are well marked and of the proper sign, although 
often with considerable scatter; while for intermediate regions the 
progression with intensity is less rapid. Between the null points 331° 
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and 61° the displacements are preponderantly positive, while be- 
tween 61° and 151° they are negative. The general conclusion is that 
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Fic. 6.—Relationship between residual radial velocity and intensity of the D lines 
for individual stars grouped according to galactic longitude. An arc above a plotted 
point indicates that the star is 10° or more north of the galactic equator; below, 10° or 
more south. The dashed line on each diagram marks the locus at which displacement 
equals one-half the equivalent width. 

Region 168°-177°.—The divergent point at the lower right corresponds to HD 46300, 
cAo. The data appear reliable. 

Region 182°-212°.—The highest point corresponds to HD 50896, Ob; the data have | 
medium weight. The point next below corresponds to HD 52382, B2; the data have 
medium weight, but H and K yield a velocity lower by 9 km/sec. 


the influence of galactic rotation is plainly manifest in these inten- 
sity-displacement diagrams. 
For lines whose chief cause of widening is galactic rotation, the 
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measured displacement should equal one-half the actual width and 
should equal or exceed one-half the equivalent width, depending on 
whether or not complete saturation occurs. In no event should the 
displacement be less than one-half the equivalent width. 

For convenient reference the locus of points for which the dis- 
placement equals one-half the equivalent width is indicated on each 
diagram by a dashed line. On only one diagram, longitudes from 
g° to 33°, does this line provide a passable representation of the 
plotted points. On the others, nearly all the points correspond to 
displacements numerically less than half the equivalent widths. This 
fact confirms the previous finding that some phenomenon other than 
rotation must play an important part in producing the observed in- 
tensities. Quite similar was Beals’s conclusions ‘‘that in many cases 
some mechanism other than galactic rotation exercises an important 
influence in the formation of the line profile.” 

A few displacements numerically greater than half the equivalent 
widths occur, notably in the region 182°-212°, where the displace- 
ments exhibit a large scatter although the range of intensity is small. 
The data appear to indicate considerable irregularity in the motions 
of interstellar sodium vapor. 
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ANALYSIS OF THE INTENSITIES OF THE 
INTERSTELLAR D LINES* 


O. C. WILSON AND PAUL W. MERRILL 


ABSTRACT 


Previous investigations of interstellar lines are briefly reviewed. 

The interstellar D-line intensities measured at Mount Wilson in the spectra of over 
200 early-type stars are critically discussed with the object of deducing the distribution 
and motions of the sodium atoms in space. If the atoms are assumed to have no rela- 
tive motions other than those due to temperature, the D-line intensities should follow 
a normal curve of growth. This assumption is shown to be incorrect because (1) the 
curve derived from the measured intensities does not have the proper shape and (2) it 
leads to very improbable distances for the stars with strong lines as compared to those 
with weak lines. 

The introduction of galactic rotation gives a good representation of the mean de- 
crease of the doublet ratio D2/D1 from 2.0 to 1.2 with increasing intensity. A more de- 
tailed examination proves conclusively, however, that the rotational velocity gradient 
cannot be responsible for the widths of the lines because, for a given distance, the line 
intensities and the doublet ratio are essentially the same for all galactic longitudes. 

The observations can be accounted for by making the hypothesis that interstellar 
sodium occurs in discrete aggregations or clouds, which, while participating in the gen- 
eral galactic rotation, have in addition considerable random motions. This assumption 
requires no dependence upon the galactic longitude and succeeds in explaining how the 
D lines can strengthen almost linearly with the distance, while the ratio of D2 to D1 is 
only 1.2. By fitting computed curves to the observational data, the following numerical 
results are obtained: an upper limit to the temperature is 44,000°; the density of sodium 
atoms in the ground state is approximately 10-3" gm/cc; the linear dimensions of the 
clouds, in the galactic plane, are of the order of 700 parsecs; and the spread of the ran- 
dom cloud velocities need not exceed 15 or 20 km/sec. 


SUMMARY OF PREVIOUS INVESTIGATIONS 


Preceding work on interstellar line intensities may, for the pur- 
pose of a brief review, be conveniently classified as either theoretical 
or observational. The exact historical order is partially disregarded 
in the following summary, which begins with the theoretical aspect 
of the problem. 

The first analysis along modern lines of the conditions prevailing 
in interstellar space was made by Eddington in his 1926 Bakerian 
lecture.t In this work Eddington used rather general physical argu- 
ments to arrive at estimates of the density and temperature of inter- 
stellar calcium and sodium, and it is noteworthy that his results did 
not depend upon accurate knowledge of the intensities of the inter- 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 


ton, No. 570. 
t Proc. Roy Soc., A, 111, 424, 1920. 
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stellar lines of these elements. Measurements of intensities had not, 
in fact, been made at the time; and only very rough estimates were 
available as an approximate check on the validity of the computa- 
tions. 

Gerasimovi¢ and Struve? in 1929 utilized eye estimates of inter- 
stellar line intensities by Struve,} in conjunction with probable 
stellar distances, to derive a reasonable value for the density of ion- 
ized calcium in space. Their determination was based upon the as- 
sumption that the strength of a line depends only on the atomic ab- 
sorption coefficient and the number of atoms in the line of sight. 

The fact that the intensity of an absorption line is a function of 
the relative atomic velocities, as well as of the atomic constants and 
total number of atoms, was made clear by the work of Schiitz in 
1930.45 Unséld, Struve, and Elvey° first applied these ideas to the 
problem of interstellar lines. They recognized that the theory of 
galactic rotation, as developed by Oort and others, would provide a 
linear velocity gradient along all lines of sight in the galactic plane 
except those pointing toward and away from the galactic center and 
at right angles thereto. If the interstellar matter participated in the 
rotation, the resultant differential velocities should be taken into 
account in computing the intensities of interstellar absorption lines. 
‘Investigation showed that the doublet ratio, K/H, could theoreti- 
cally have any value in the range from 1 to 2, depending upon the 
distance and position of the star in whose spectrum the lines oc- 
curred. The amount of observational material included in the in- 
vestigation was insufficient to provide a satisfactory comparison 
with the theoretical results. 

More recently Eddington’ also has adopted the idea that galactic 
rotation may be chiefly responsible for the widths and total absorp- 
tions of interstellar lines. He combined a linear velocity gradient 
with the usual atomic absorption coefficient and computed the in- 
tensities of the resultant lines. His analysis, which treats the density 
of the interstellar gas as a variable parameter, indicates that after 
the density has attained a value sufficiently great to produce a 


2 Ap. J., 69, 7, 1920. 5 Struve and Elvey, Ap. J., 79, 409, 1934. 
3 Ibid., 67, 353, 1928. 6 Zs. f. Ap., 1, 314, 1930. 
4Zs.f. Ap., 1, 300, 1930. 7M.N., 95, 2, 1934. 
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saturated line, the density may then be increased about five hundred 
fo'd before the line shows a further increase in total absorption. 
Eddington considers this result to be the explanation of the hitherto 
anomalous fact that the Nai and Ca 11 lines are observed to be of 
approximately equal strength, while his previous purely theoretical 
calculations had shown that the ionized calcium atoms in space 
should be about three hundred times as numerous as those of neutral 
sodium. By considerations of this nature and by comparison with 
measured line intensities, Eddington was led to estimates of the 
densities of interstellar calcium and sodium which, owing partly to 
insufficient observational data, cannot be considered as completely 
satisfactory. Among other results of this stimulating paper, there is 
one which is very easy to test observationally. This is the fairly 
obvious statement that, if for a saturated line the width and dis- 
placement are both due to a linear velocity gradient along the line 
of sight, the radial velocity cannot exceed half the equivalent width 
expressed in velocity units. 

Struve was the first to make a systematic attack on the problem 
from the observational side. He made eye estimates of the inten- 
sities of the interstellar Ca 11 lines in the spectra of a large number of 
early-type stars. Although his individual values were undoubtedly 
of low weight, owing to the inhomogeneity of his material and the 
inherent uncertainties of such work, they were of great statistical 
importance. He found that the intensities of the lines increased with 
the distances of the stars in whose spectra they were observed. This 
fact has been amply confirmed by later work, although there is evi- 
dence that the precise correlation between intensity and distance 
found by Struve is probably not correct. 

The extensive and well-known researches of Plaskett and Pearce*® 
had as their main objective the determination of the constants of 
galactic rotation from measures of the displacements of the stellar 
and interstellar lines in the spectra of many early-type stars. These 
authors made careful eye estimates of the intensities of numerous 
interstellar Ca 11 lines, by comparison with a set of standard artificial 
lines, the chief purpose being to arrange the stars into groups on the 


8 Pub. Dominion Ap. Obs., 5, 167, 1933. 
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basis of the strength of the K line. The intensities derived in this 
manner are probably representative of the highest accuracy that 
has been attained by visual comparison. Nevertheless, the superi- 
ority of photometric measurement of line strength over the best 
previous results is clearly brought out by E. G. Williams’ comparison 
of his own measures with the Victoria values.’ He found a consider- 
ably larger scatter on plotting the two sets of data against each other 
than would be anticipated from the size of his probable errors alone. 
Evidently the greatly increased labor of measuring interstellar line 
intensities by the usual photometric methods is entirely justified, 
particularly where accurate determinations of the relative strengths 
of the two members of a doublet are required. 

Williams was concerned chiefly with determining the relationship 
of interstellar K intensity to color excess; and his results do not, for 
the most part, have a direct bearing on the present investigation. It 
is noteworthy, however, that he encountered considerable difficulty 
in obtaining reliable values of the doublet intensity ratio, K/H, 
owing to the blending of the interstellar H line with stellar He. We 
have seen that accurate knowledge of this ratio is very desirable for 
comparison with theoretical predictions, and the much greater free- 
dom of the sodium D lines from blends with stellar lines is a strong 
argument in favor of employing them in an extensive investigation. 

In his most recent paper Beals” has collected and discussed his 
published observational material on both calcium and sodium and 
has arrived at certain important results. He finds, on plotting radial 
velocities of interstellar lines against their intensities, only a weak 
correlation at best, and concludes that Eddington’s simple picture 
of line formation cannot be adequate. He suggests irregular group 
motions in the interstellar cloud as the main widening agent, at least 
for the near-by stars, and in support of this idea finds that high-dis- 
persion spectrograms of three bright B-type stars show the K lines 
to consist of two components, neither of which, judged by velocity 
or appearance, is at all likely to be of stellar origin. He investigates 
also the doublet ratios, K/H and D2/D1, and finds good evidence 
that in both cases the ratio attains the value 2 for weak lines and de- 


9 Mt. W. Contr., No. 487; Ap. J., 79, 280, 1934. 1 M.N., 96, 661, 1936. 
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creases as the lines become stronger. His data are probably not com- 
plete enough to determine accurately the variation of the ratio with 
line strength, since there appears to be considerable scatter. 


DISCUSSION OF MOUNT WILSON DATA 
The need for further measures of interstellar line intensities has 
been made clear by the foregoing summary of previous work in this 
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Fic. 1.—Relationship between equivalent widths, D1 and D2 


field, and the advantages to be derived from an extensive study of 
the sodium D lines have been pointed out. We have therefore under- 
taken the measurement of the intensities of the interstellar D lines in 
the spectra of over 200 stars, well distributed over the range of 
galactic longitude observable from Mount Wilson. Our measures, 
together with such relevant matters as discussions of the photo- 
metric technique and general accuracy, will be presented fully in 
another Contribution." In the present paper the data are analyzed 


™ Mt. W. Contr., No. 576; in press. 
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in an attempt to deduce the temperature and density of sodium in 
interstellar space and the mechanism of interstellar line formation. 
The analysis consists of a comparison of the observational data with 
the results of computations based upon various assumptions as to 
the motions and distribution of the interstellar gas, the simplest 
hypothesis being considered first, and more complicated ones in- 
troduced as the need for them arises. 

In Figure 1 the total absorption of D1 is plotted against that of 
D2 for each star.” This diagram illustrates the general consistency 
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Fic. 2.—The doublet ratio, D2/D1, as a function of the intensity of D2. Circles 
correspond to stars near the null directions of galactic rotation. 


of the measures and shows clearly that the ratio of D2 to D1 tends 
to be decidedly larger for the weaker lines than for the stronger ones. 
This interesting feature of the data is best exhibited by plotting the 
doublet ratio, D2/D1, against D2, as has been done in Figure 2. 
The curves drawn in Figure 2 will be referred to in a later section; 
for the present, attention is directed to the distribution of the plotted 
points. It is evident that, for the weakest lines, the doublet ratio is 

" Hereafter the symbols D1 and D2 will frequently be used to signify the total ab- 


sorptions of the lines as well as to refer to them in the more general sense. The meaning 
in any case will be clear from the context. 
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about 2; and, as the lines increase in strength, the ratio rapidly di- 
minishes to about 1.2 and remains constant throughout the range of 
intensity covered by the observations. Beals’s data are in essential 
agreement with these results. 

The correlation of high ratio with low intensity suggests the pos- 
sible blending of interstellar D2 with some other weak line. This 
point was forcibly brought to our attention by the fact that several 
stars with weak D lines were found to yield doublet ratios definitely 
in excess of the theoretical upper limit, 2. The anomalous ratios in 
these spectra were traced to a blending with D2 of the stellar C 
pair at 5889.97 and X 5891.65. The presence of these C 1 lines 
can apparently be inferred from the behavior of the other well- 
known C 11 lines at \ 4267 and A 3919. If the latter lines can be seen 
easily, the doublet ratio should probably not be used unless the 
equivalent width of D2 exceeds 0.15 A. If D2 is greater than o.15 A, 
the contribution of the C 1 lines, if present, will probably not be 
sufficient to produce a serious distortion of the doublet ratio. This 
criterion has eliminated several stars from consideration. Figure 2 
includes only those stars in which we feel that the ratio is reasonably 
free from the influence of the stellar C 11 lines." 

The appearance of Figure 2 suggests an attempt to explain the 
doublet ratios by means of a normal curve of growth. The first step 
is to construct a curve of growth from the measured intensities. If we 
assume a homogeneous distribution of atomic velocities throughout 
space, the construction of the curve from the measured intensities 
would be a very simple matter provided the number of atoms in the 
line of sight to each star were known. But knowledge of the total 
number of atoms to any star requires knowing its distance, as well 
as the density of interstellar sodium; and, since we hope to deter- 
mine these data from our measures, we do not wish to assume them. 

First consider the physical conditions under which a normal curve 
of growth would be anticipated. These are typified by a laboratory 
experiment in which a tube contains atoms whose only motions are 
those due to the temperature, 7, assumed constant throughout the 


13 We are indebted to Dr. J. A. Pearce, of the Dominion Astrophysical Observatory 
for data on the violet lines of C 1 in a number of spectra. 
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tube. If, in such an experiment, the number of effective atoms, Nf, 
can be varied and the corresponding equivalent widths, W, meas- 
ured, then the curve resulting when log Nf is plotted against log W 
is the curve of growth for the given line at the specified value of T. 
Such curves are characterized by the fact that they represent rela- 
tions of the form W « (Nf)*, where the exponent a can range from 1 
to nearly o. Curves of growth computed for sodium are shown by 
the dashed lines in Figure 3. The straight line on the right represents 
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Fic. 3.—Curves of growth for the D lines. The continuous curve is constructed from 
the measured equivalent widths. The dashed curves are theoretical. 


the case where the atoms are quiescent and the absorption lines are 
produced by pure radiation damping. This line corresponds to 
a = 3. In general, for a given 7, the EW’s will start on the left- 
hand line, where a = 1, rise to a point determined by 7, and then 
will follow a curve such as those shown, ultimately joining the 
square-root branch on the right. The important case where W is not 
a function of Vf alone, but of the distance as well, will be discussed 
later. 

Suppose now that we have measures of the total absorptions of a 
doublet, such as the D lines, for a series of values of the number of 
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atoms, V. Let the successive values of NV be so close together that, 
if Wi and Wi are the EW’s of the ith pair, we have 


<wWi. 


Any pair of values corresponds to effective numbers of atoms in the 
ratio 1:2, since this is the ratio of the f values for the D lines. Hence, 
if log W; and log W? are plotted against log N and log N + Alog N, 
respectively, where JN is arbitrary and A log N = log 2 = 0.30, the 
straight line joining these two points will have a slope which is equal 
to that of the curve of growth at some point between W; and W3. 
Successive values of log W, and log W, are then plotted against the 
same values of log NV, and each pair of points is joined by a straight 
line. Thus we have a series of line segments of varying inclination, 
lying one above the other and each having a slope equal to that of the 
curve of growth at some point. Next, the second of these segments 
is moved horizontally to the right until it touches the first; then the 
third until it touches the second; etc. Finally, if a smooth curve is 
drawn such that all the segments are tangent to it, we should have 
a good approximation to the true curve of growth. This method was 
tested by reading pairs of points from a computed curve of growth 
and using them in the foregoing manner to reconstruct the original 
curve. The agreement between the two curves was quite good, and 
the method outlined appears satisfactory. 

This scheme was applied to the measured EW’s of interstellar 
lines. All the stars were arranged in order of increasing D1, and 
means of D1 and D2 were taken in groups of ten, overlapping five 
each time. These mean values of D1 and D2 satisfied the condition 
noted above and were used to construct the curve of growth for the 
interstellar sodium lines. Finally, this curve was superposed on the 
theoretical ones and shifted horizontally to obtain the most generally 
satisfactory fit. The result is shown by the solid curve in Figure 3. 
Evidently the observed curve corresponds to a value of the param- 


eter 6 in the neighborhood of o.1 or 0.2, where, in the notation of 


Struve and Elvey,' 6 is defined by 6 = \,v./c, with v, equal to the 
most probable velocity of the atoms. 

By fitting the observed curve to the theoretical curves, we have 
supplied it with a true log NV scale. Thus, if the lines were formed 
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under conditions which yield a normal curve of growth, Figure 3 
would give the number of atoms active in producing any given ob- 
served line. The assumption of uniform sodium density then leads 
immediately to the relative distances of the stars. Assuming the 
legitimacy of this process, we note from the figure that there are 
about three hundred times as many atoms involved in the produc- 
tion of our strongest lines as in the weakest. Since there is good rea- 
son to believe that the weakest lines correspond to stellar distances 
of at least 100 parsecs, the stars with intense lines should be 30,000 
parsecs distant, whereas their actual distances are probably only 
two or three thousand parsecs. This difficulty could be avoided, of 
course, by postulating nonuniform density; but, when it is considered 
in conjunction with the fact that the observed and theoretical curves 
definitely disagree in shape, the proper conclusion appears to be that 
interstellar line intensities are not represented by a normal curve 
of growth. Accordingly, the solution shown in Figure 3 is rejected. 

Since the hypothesis of uniform density, together with that of 
pure temperature motion or small-scale turbulence, appears inca- 
pable of giving a satisfactory solution, motions dependent on the dis- 
tance, as suggested by the theory of galactic rotation, may next be 
considered. The appearance of Figure 2 suggests the following ten- 
tative ideas. For the near-by stars the temperature or other small- 
scale motions of the atoms are the determining factor in producing 
the absorption lines, while for greater distances the effect of rotation 
will predominate. Hence the proper manner of introducing rotation, 
without sacrificing the possibility of deducing both the temperature 
and density of the sodium, appears to be the following: 

Consider an element di of any line of sight having a cross-section 
1 sqcm. The optical thickness‘of this element for wave-length X is 
given 
dr, = an, (1) 

where e, m, c, and \, have their usual meanings, ” is the number of 
active atoms per cubic centimeter, f the oscillator strength, and 
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14 A. Unséld, Mt. W. Contr., No. 377; Ap. J., 69, 209, 1929. 
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on the assumption of a Maxwellian velocity distribution correspond- 
ing to temperature 7. We assume and 6 constant along the line 
of sight. If the radial velocity of the element is taken to be propor- 
tional to its distance, we may write 


dr, = a = dr, , 


whence, 


where X, and X, are the values of \, at the two ends of the line of 
sight. Setting (A, — A)/b = x, we have 


Finally, convenient units are introduced, and the expression for 7 
becomes 


.98 , AX 


4 x40 VI (2) 


In equation (2) the various symbols mean: 


n = number of active atoms per cubic centimeter. 
V = velocity gradient along the line of sight in kilometers per second per kilo- 
parsec. 


dete _ [2kT 


1 = length of light-path in parsecs. 
AX = wave-length distance in angstroms to the line center. 
¢ = probability integral. 


This expression for rt is quite complicated and does not yield a 
simple result for the EW as does that derived by Eddington, but it 
has the advantage of taking accurate account of temperature mo- 
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tions as well as galactic rotation. To apply equation (1) the EW is 
given, as usual, by 


w= (3) 


with 


(4) 


The method of computing line intensities from equations (2), (3), 
and (4) deserves some comment. Originally the contours derived 
from (2) and (4) were plotted and integrated with a planimeter. This 
is a time-consuming process when the number of lines is fairly great, 
and a much shorter method was found to be the following: Let r, be 
the value of r at the center of the line, i.e., for AA = o. Then for any 
line we may always write 


W = 2(1 — 7r-)AXm 


where A),, is the half-width at some point within the line. For the 
kind of lines given by (2) and (4) it turns out that the correct value 
of AX,, isthe AX at r = (1 + 7.)/2, i.e., at the point halfway from the 
continuous background to the bottom of the line. Values of W ob- 
tained in this manner have been compared with those derived by 
planimeter measures over the entire range of intensity needed for 
this investigation and always agree within a per cent or two. In- 
cidentally, this method is applicable a!so to the case where 7 is given 
by (1) and is shorter in practice than the series expansions used by 
Struve and Elvey.s 

There now remains the step of comparing the observed and com- 
puted line intensities in such a way as to determine as many as pos- 
sible of the unknown parameters entering into equation (2). This 
comparison might be effected in various ways. We have adopted the 
following procedure for reasons which will presently become clear. 

A little thought will show that, unless additional data are given, 
measures of line intensity alone can scarcely determine the velocity 
gradient V. This quantity is most logically found from velocity meas- 
ures, and we have therefore adopted the value V = 17 km/sec/kilo- 
parsec in the direction of a maximum point of the rotational effect. 


D2/D1 
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This is the value assumed by Eddington in his discussion and is 
very close to the final result of Plaskett and Pearce. In order to use 
all our stars statistically, we must take the appropriate mean value 
of the gradient, which is 2/7 times that at a maximum, or 10.8 
km/sec/kiloparsec. V is now set equal to 10.8 in equation (2), and 
the total absorptions of D1 and D2 are computed for a set of stand- 
ard distances ranging from 100 to 6000 parsecs and for various com- 
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Fic. 4.—(a) Illustration of the method of determining and 6 by comparing the 
results of computation with the observations. Dots are observed normal points. For 
the continuous curves, b=0.05; for the dashed curves, 0.15. The upper curve of each 
pair is for m=0.95X 107°; the lower one for n=3.8X10°°. 

(b) Representation of observed intensities (and ratios) by the accepted values, 
n=2.0X10% atoms per cubic centimeter, b=o.11 A. Distances in parsecs are in- 


dicated on the curves. 


binations of the two parameters m and b. From these computed 
values a series of curves of D2/D1 against D1 are plotted, and m and 
b are adjusted for the best possible fit with the corresponding ob- 
served curve. This method proves to be more sensitive in practice 
than might be judged from the foregoing description. 

The process is illustrated in Figure 4a. The dots are normal 
points representing the observed relation for the sodium lines. The 
two solid curves are for 6 = 0.05, and the broken ones for 6 = 0.15. 
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The uppermost curve of each pair corresponds to n = 0.95 X 107°, 
and the lower one to m = 3.8 X 10-’ atoms per cubic centimeter. 
The fit of the observed with the computed curves depends chiefly 
upon » for large values of D1 and chiefly upon 6 for small values, 
thus permitting semi-independent estimates to be made of both 
quantities. Further computation resulted in the curves of Figure 40. 
The values of the parameters entering here are m = 2.0 X 10~® and 
b = 0.11. The three curves are, from top to bottom, for values of 
V equal, respectively, to 17.0, 10.8, and o km/sec/kiloparsec, cor- 
responding to conditions at a maximum point, an average point, 
and a null point. The observations should fit the upper curve some- 
what better than the middle one for the stronger lines, since it hap- 
pens that we have observed more of these stars near maximum points 
than near null points. So far, the results of this analysis seem satis- 
factory, and the derived numerical values of the constants are rea- 
sonable. The temperature corresponding to b = 0.11 A is 44,000°, 
rather higher than is usually assumed; but this value may, of course, 
include some small-scale turbulence, in which case the true tempera- 
ture would be less than the foregoing value. The density of neutral 
sodium comes out to be 7.6 X 107%? gm/cc, in quite close agreement 
with Eddington’s estimate. 

We have thought it of sufficient general interest to carry through 
these computations based on the galactic rotation hypothesis, al- 
though a simpler approach to the problem’ *’ has already thrown 
considerable doubt on the validity of the assumption that rotation 
is the main cause of the widening of interstellar lines. Thus far, only 
the statistical aspect of the problem has been considered, and we 
must now investigate in detail whether or not the dependence of 
line intensity on galactic longitude is in agreement with the predic- 
tions of rotation theory. 

This dependence on galactic longitude could be examined in a 
variety of ways. The method chosen is based on the desire to use 
the derived constants, ” and b, for the determination of more accu- 
rate distances of individual stars from measures of intensity of in- 
terstellar sodium lines. Evidently it is now possible to compute from 


15 Mt. Wilson Contr., No. 569; Ap. J., 86, 28, 1937. | 
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equations (2) to (4) a series of curves relating intensity and dis- 
tance at various galactic longitudes. Distances obtained in this way 
may then be compared with those derived by the ‘‘spectroscopic”’ 
method. The latter values may, of course, involve systematic errors; 
but such errors should not be functions of galactic longitude. Hence, 
if these two sets of distances are found to bear different relations to 
each other for different longitudes, the fault must lie, almost cer- 
tainly, with equation (2) and the underlying theory. 

One set of distances was computed in the usual way from the ap- 
parent magnitudes and assumed absolute magnitudes.® Another set 
was determined by the measured line strengths, and the curves com- 
puted by means of equations (2) to (4). The stars were divided 
into three groups for which A/, the longitude difference between the 
star and the nearest null point, lay between o° and 11°, 12° and 23°, 
24° and 45°, respectively. The longitude of the galactic center was 
taken to be 331°. The two values of the distance were plotted against 
each other for every star, with the results shown in Figure 5a, 8, c, 
where abscissae are distances based on line intensity and rotation 
theory, and ordinates are distances derived from the spectral types 
and apparent magnitudes. Obviously, the relation between the two 
sets of distances is very decidedly a function of the galactic longitude 
of the stars considered. For those stars not near rotational null 
points, the distances computed from observed line intensity are in 
reasonable agreement with the other set. On the other hand, for the 
stars nearest the nulls, where the rotational velocity gradient is 
small, much greater distances must be assumed in order to yield in- 
tensities equal to those observed. These facts may be put into other 
words, as follows: If galactic rotation plays the chief role in the 
formation of the interstellar sodium lines, then, near the null points, 
where the rotational effect is small, the lines should be weaker, for a 
given distance, than near the maxima. Observationally this is not 
found to be true, for the lines have much the same strength, for a 
given distance, at all longitudes. Clearly, some potent widening 
effect other than rotation must be sought. 

On the entirely groundless assumption that the rotational velocity 

© For details see Mt. W. Conir., No. 569; Ap. J., 86, 28, 1937. Especial care was 
taken to avoid errors in spectral type. 
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gradient V = 17 km/sec/kiloparsec applies to stars at all longitudes, 
we find the result shown in Figure 5d, which, while having no signifi- 


Parsecs 
AL at 
To 1° b 12 TO 23° 
5000 |- 
4000 }- 
3000 }- 3 
~ 
=| 
a 
ee 8 e ? e 
o 
at 
24° To 45° d To 45° 
.& 5000 - 
= 
Q 
n 
3000 a e oe e ee 
1000 
° 2000 4000 ° 2000 4000 parsecs 


Distances from spectral types 


Fic. 5.—Comparison of stellar distances derived from measured D-line intensity 
and rotation theory (ordinates) with those deduced from spectral types and apparent 
magnitude (abscissae). A/ is the longitude difference, regardless of sign, between a star 
and the nearest null point of galactic rotation. In the lower right-hand diagram the 
velocity gradient, 17 km/sec/kiloparsec, is assumed to apply to all stars regardless of 


galactic longitude. 


cance on the basis of rotation theory, is likely to be of considerable 
importance in determining the velocity distribution of the inter- 


stellar sodium. 
Thus we find, as did Beals, that rotation is inadequate as the 
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source of the strength of interstellar lines. However, one more set 
of diagrams proves interesting in this connection. In Figure 6 are 
plotted residual interstellar radial velocities against total absorp- 
tions in kilometers per second, the grouping of stars being the same 
as that of Figure 5.'7 Evidently, for those stars nearest the null 
points a correlation between these quantities could be found only 
with the exercise of considerable imagination. For the intermediate 
group the situation is not much better, except possibly for a half- 
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Fic. 6.—Residual radial velocities from the D lines compared with the mean equiva- 
lent widths in km/sec, for groups of stars used in Fig. 5. 


dozen or so of the stars with strongest lines. But for the stars in the 
vicinity of the maxima of rotational effect, a correlation is un- 
doubtedly present. In all three diagrams a very great scatter is evi- 
dent. Lines having two or three times their maximum allowable 
strength, on the basis of rotation, as judged by their measured veloc- 
ities, are numerous. The measures of both velocity and intensity are 
sufficiently accurate to insure that this scatter is real for the most 
part. 

At this point it is well to bring together the facts we have ascer- 
tained and to indicate the requirements of a successful theory. We 
have found: 

17 The data are plotted in greater detail in Mt. W. Contr., No. 569, Figs. 5 and 6; 
Ap. J., 86, 41, 42, 1937. 
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1. The doublet ratio D2/D1 is approximately 2 for our weakest 
lines, thereafter decreasing and maintaining the value 1.2 for the 
stronger lines. 

2. The attempt to explain the observations on the basis of homo- 
geneous conditions in interstellar space is considered unsuccessful 
because (a) the shape of the observational curve differs from the 
theoretical by a larger amount than we feel is to be attributed to ob- 
servational error, and, (b) on this hypothesis many stars must be 
placed at improbably great distances. 

3. By introducing galactic rotation a good quantitative explana- 
tion of the variation in the doublet ratio is obtained, and reasonable 
values for the temperature and density of the interstellar sodium can 
be deduced by using the measured intensities statistically. 

4. The simple rotational picture becomes untenable, however, 
when the behavior of line strength as a function of galactic longitude 
is investigated. For a given distance, the lines are not markedly 
weaker near the null points than near the maxima. Moreover, the 
ratio D2/D1 shows no dependence on galactic longitude. This fact 
is illustrated in Figure 2, where the stars are divided into two groups, 
depending on whether they lie closer to a null point or to a maximum 
point. 

5. The presence of the rotational effect is, however, indicated by 
the fact that the correlation between velocity and intensity is much 
better for stars near the maxima than for those near the nulls 
(Fig. 6). Apart from this, the great scatter in Figure 6 is evidence 
that some powerful widening agent besides rotation is at work. 

A further attack on the problem, which seems to have proved suc- 
cessful, will now be presented. As a preliminary step, attention is 
directed to some of the diagrams in the regional study of the sodium 
lines by one of the authors.’ Figure 1 of that paper shows, for the 
stars in certain small areas of the sky, the mean D-line intensity 
plotted against distance. The distances used in these figures are 
those derived from the spectral types and apparent magnitudes. In 
practically all cases the intensities increase with the distances; but, 
while the rate of increase is definitely not the same for all regions, 
the scatter of the points for a given area from the mean relation for 
that area is generally fairly small. Combination of the results for 
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twenty-three areas leads to Figure 2 of the regional paper, in which 
the mean curve has been corrected for an assumed space absorption 
of 0.35 mag. per kiloparsec. The figure shows that, in the mean, the 
intensities of the interstellar sodium lines increase very nearly pro- 
portionally to the distances of the stars in whose spectra they are 
observed. 

The mean sodium-line intensities are plotted against the corrected 
distances for all stars within 10° of the null points in Figure 7 of the 
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Fic. 7.—Mean D-line intensity plotted against distance for stars within 10° of the 
rotational null points. Crosses represent stars near longitude 331°, circles those near 
151°, and dots those near 61°. The curves are explained in the text. 


present paper. Although effects of galactic rotation should be very 
small, the intensities nevertheless are again nearly proportional to 
the distances. The stars at the three observed null points have been 
given distinguishing symbols, and it appears that the relationship of 
intensity to distance is not quite the same at the three nulls. The 
stars near galactic longitude 331°, for instance, seem to lie somewhat 
lower in the diagram than those near longitude 61°. Such a difference 
could be accounted for by postulating somewhat heavier local ob- 
scuration at 331° than at 61°. The group of stars with strong lines 
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at distances less than 800 parsecs is somewhat more troublesome. 
Two of them, a Cygni and v Cephei, indicated in the figure, are of 
type cA2, and there is evidence that such stars are intrinsically 
brighter than was assumed in assigning them distances. Hence these 
two stars are perhaps plotted too far to the left. The point between 
them corresponds to HD 203025, a nebulous star whose sodium lines 
so greatly exceed the K line in strength as to indicate almost cer- 
tainly that the D lines are abnormally strong. The remaining stars 
in this group seem to be normal. 

These sources of error appear, however, not to explain all the 
scatter in Figure 7. In order to bring a number of the plotted points 
into agreement with the mean, it must be assumed that either the 
intensities or the distances, or both, are incorrect by rather large 
factors, of the order of 2 in some instances. The authors believe it 
to be almost certain that a considerable part of this scatter is real, 
since the analogous diagrams for small regions previously referred to 
indicate the purely observational scatter to be considerably less. 

We may now set forth the facts which must be met by any suc- 
cessful explanation. These are: 

1. The decrease in doublet ratio from 2.0 for weak lines to 1.2 for 
stronger lines (Fig. 2). 

2. The almost linear increase of intensity with distance in the null 
directions without accompanying bodily displacements of the lines."® 

3. The scatter of the observations as shown in Figure 7. 

4. The observed displacements of the lines, which in general are 
in good agreement with the predictions of the rotation hypothesis. 

If it is agreed that the widening and displacement of interstellar 
lines are manifestations of Doppler effect, the only reasonable hy- 
pothesis satisfying the requirements is that the interstellar gas exists 
in discrete masses, or clouds, and that these clouds, while partaking 
of the general galactic rotation, have additional relative motions. 

The usefulness of the notion of discrete clouds having a spread of 
velocities is made apparent by considerations like the following. 
Suppose that a star is observed at such a distance that its light has 
traversed one cloud. The interstellar sodium lines in its spectrum 


8 The displacements are plotted against longitude by Plaskett and Pearce; see foot- 
note 8. Measurements of interstellar sodium lines yield essentially similar results. 
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will have intensities, D1 and D2, and a ratio, D2/D1, which will 
depend upon the dimensions of the cloud and the density and tem- 
perature of the sodium within it. Now the light of a second star at 
about twice the distance of the first may be considered to pass 
through two clouds, whose radial velocities differ sufficiently to pre- 
vent blending of the two sets of lines."® Consequently, the inter- 
stellar lines in the spectrum of the second star will have about twice 
the intensity of those in the first, but the doublet ratio will be the 
same for both stars. As still more distant stars are observed, a point 
will be reached ultimately where overlapping of the individual cloud 
lines becomes prevalent, and from there on the intensities increase 
less rapidly, approximating the behavior of the flat portion of a nor- 
mal curve of growth having a parameter P determined by the aver- 
age velocity of the clouds. 

The preceding example should help to clarify the following more 
general discussion. The cloud hypothesis makes four parameters 
available for adjustment. These are: n, the average number of 
ground-state sodium atoms per cubic centimeter within the clouds; 
b, the “temperature” parameter expressing the average kinetic mo- 
tions of the atoms; P, the corresponding quantity expressing the 
spread of the cloud velocities; /,, the average length of a line of sight 
within a cloud. The problem is to choose such values of these quan- 
tities that the line intensities computed from them most nearly agree 
in behavior with the observed values. In other words, the observa- 
tional data shown by the plotted points in Figures 2 and 7 must be 
fitted as accurately as possible by computed curves, and the achieve- 
ment of a satisfactory fit fixes the values of the parameters. The 
process is thus quite analogous to that already employed in this 
paper in the attempt to satisfy the observations by simple rotation 
theory, and is, essentially, a method of trial and error. 

The authors are greatly indebted to Dr. R. Minkowski for sug- 
gesting and carrying out the computations necessary for a compari- 


19 Tt turns out in practice that a difference in wave-length of about o.1 A, corre- 
sponding to a velocity difference between the clouds of about 5 km/sec, is sufficient to 
prevent serious overlapping of the lines. Such a small difference is, however, not enough 
to separate the lines with the relatively small dispersion used in this investigation, and 
consequently the interstellar lines appear single on our spectrograms. 
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son of the hypothesis with the observations. His solution, which was 
greatly facilitated by the use of tables published by Ladenburg and 
Levy,” is exhibited by the curves shown in Figures 2 and 7. 

The first step is to reproduce that portion of the ratio-curve in 
which the doublet ratio decreases from 2 to 1.2, and this amounts 
to choosing values of b and of the average total number of atoms 
traversed by the line of sight in the cloud nearest the observer. In 
Figure 2, the solid curve between points A and B, which has been 
taken from a curve of growth with b = 0.11, is considered to be in 
satisfactory agreement with the observations. Turning now to Fig- 
ure 7, we find that the same curve of growth gives the solid curve 
between A and B for the relationship between mean D intensity and 
distance. The point B is seen to lie at a distance of 720 parsecs. 
Hence, dividing the average total number of atoms in the line of 
sight in the first cloud by the number of centimeters in 720 parsecs 
yields the value of . 

After passing points B in the diagrams, the average line of sight 
is assumed to enter a second cloud having the same physical char- 
acteristics as the first one, but differing sufficiently from the latter 
in radial velocity to prevent overlapping of the two sets of absorp- 
tion lines. Thus, adding the intensities of the second set of lines 
carries the two curves from points B to points C. Similarly, a third 
cloud takes the curves to points D, beyond which a scarcity of ob- 
servations renders further calculation futile. 

The argument as presented above corresponds to the idealized 
case where the clouds, all identical, are penetrated centrally one 
after the other by the line of sight. A little thought will show that, 
when the more likely case of only partial penetration of some of the 
clouds is considered, the theory demands a considerable scatter in 
diagrams like Figures 2 and 7. Suppose, for instance, that we are 
somewhere within a cloud and that in some direction we find its 
boundary at a distance of only 250 parsecs, corresponding to the 
point designated £ in the figures. After passing this point, the sec- 
ond cloud is entered, and the variation of intensity with distance will 
follow the dashed curves from E to F. Similarly, gaps between the 


20 Zs. f. Phys., 65, 189, 1930. 
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clouds will produce scatter. If, for instance, there were a space of 
200 parsecs between the first and second clouds, the curve should 
pass through the points ABB’C’ in Figure 7 rather than as drawn. 
Thus it is clear that the present hypothesis is not only capable of 
accounting for the scatter in such a diagram but, in fact, actually re- 
quires that scatter exist. 

The numerical values of the parameters used for the calculation 
of the curves in the diagrams are of great interest. It has already 
been stated that b was taken equal to 0.11 A and that the average 
length of a line of sight through a cloud turned out to be 720 par- 
secs. The corresponding value of is 2.9 X 10-* atoms per cubic 
centimeter. The observations suffice to determine only an approxi- 
mate lower limit to the value of P, the parameter which expresses 
the distribution of cloud velocities. A glance at Figure 7 shows the 
basis for this statement. The dashed curve in the figure labeled 
P = 0.3 is computed from a curve of growth with the parameter P 
equal to 0.3 A. Clearly, as the number of clouds in the line of sight 
increases, a point will be reached where overlapping of the individual 
cloud lines will become prevalent, and from there on the rate of 
growth of intensity with distance will be slower and will, in fact, 
approximate more and more closely to that computed from a curve 
of growth having the parameter P corresponding to the distribution 
of cloud velocities. The smaller the value of P, the smaller will be 
the distance at which this effect sets in. Evidently, from Figure 7, 
P cannot be taken as much less than 0.3 without giving a disagree- 
ment with the observations. The true value of P could be found only 
by observing stars at sufficiently great distances to fix the point at 
which the rate of growth of the lines with distance begins to deviate 
appreciably from that for the nearer stars. 

There is, of course, some latitude in the values of constants de- 
termined by a method of the kind employed here, but it proves to 
be surprisingly small. In fact, it is found that must lie between 
2.9 X 107° and 3.6 X 10~°, and 6 between 0.11 and 0.09 A in order 
not to spoil the fit with either of the observed curves. It will be re- 
called that simple rotation theory yielded the values m = 2.0 X 107° 
and 6 = 0.11 A, in remarkably close agreement with the present 
ones. The two values of } should agree closely, since they are de- 
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termined in essentially the same manner; and so we adopt b = 0.11 
as the best value obtainable from our data. The values of m are 
easily reconciled since, in the rotation solution, it was assumed that 
b = 0.11 applied at all distances. It is now clear that for the more 
distant stars allowance should have been made for the distribution 
of cloud velocities. When two Maxwellian velocity distributions 


b = o.11 and P = 0.3, 8 ~o.3. If this value is used in place of 0.11 
and 7 is set equal to 2.9 X 10-°, a repetition of the rotational calcu- 
lation yields the doublet ratio 1.22 for the more distant stars, in good 
agreement with the observations. 

It remains now to discuss a few additional points concerning the 
combination of the rotational and cloud hypotheses. There are two 
distinct ways in which the clouds may conceivably participate in 
the general galactic rotation. The first case is that in which the 
whole cloud moves with the velocity of its center of gravity. In this 
event it is clear that toward a maximum, rotation would not play 
any role in producing the equivalent widths of the absorption lines, 
at least out as far as the extreme distances in Figure 7. Successive 
cloud lines would, of course, be gradually shifted in one direction, 
thus causing a displacement of the resultant line as a whole. A linear 
velocity gradient would not, however, tend to increase overlapping 
of the individual lines. On the contrary, it would, at sufficient dis- 
tances, have the effect of retarding the onset of overlapping; and 
hence, toward the maxima, a linear relationship between intensity 
and distance should continue to hold at greater distances than in the 
null directions. Our observations do not extend to sufficiently dis- 
tant stars to show this effect, if it exists. It would appear, then, that 
if the clouds move as units with the velocities of their centers of 
gravity, it was at least partly a fortuitous circumstance that caused 
the application of simple rotation theory to yield an approximately 
correct value of n. 

The other possibility is that the rotational gradient is effective 
within the individual clouds. In this case the cloud lines are widened 
and caused to overlap, and our previous solution, with the appro- 
priate value of 8, would have physical significance. Of course, if 


| 
having the parameters b and P are combined, the resultant is a 
Maxwellian distribution with the parameter 8 = Vb? + P?. With i 
i 
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differential rotation takes place within the clouds, there is the logical 
difficulty of accounting for their continued existence; but as a purely 
observational matter, the present data do not seem capable of de- 
ciding between the two possibilities. 

One important consequence of the correctness of the hypothesis 
advanced here would be that Eddington’s explanation of the equality 
of the interstellar calcium and sodium lines could no longer be con- 
sidered valid. The reason is that, to account for the doublet ratio 
for weak lines, we have found it necessary to make use of that por- 
tion of the curve of growth in which the intensity varies directly with 
the number of atoms. Hence, if the ionized calcium atoms were 
three hundred times as numerous as those of neutral sodium, the 
expectation would be that for near-by stars the H and K lines should 
be considerably stronger than the D lines. Conversely, since the ob- 
servations indicate approximate equality, the conclusion is that the 
numbers of the two kinds of atoms are roughly the same. 

To summarize briefly, we have found that the intensities and dis- 
p:acements of the interstellar sodium lines can be explained by the 
assumption that the sodium occurs in large clouds having average 
linear dimensions of the order of 700 parsecs. The clouds participate 
in galactic rotation and in addition have a velocity distribution of 
their own whose spread need not exceed a value of the order of 15 
or 20 km/sec. The temperature of the gas within the clouds cannot 
be specified exactly since any small-scale turbulence will necessarily 
be included. An upper limit to the temperature is about 44,000°, 
corresponding to 6 = 0.11 A. The best value of the average num- 
ber of neutral sodium atoms per cubic centimeter is of the order of 
3 X 107°, corresponding to a density of about 10-3" gm/cc. 


The outstanding task in the interpretation of interstellar sodium 
lines is to reconcile the low ratio of D2 to D1, 1.2 (except for very 
weak lines), with the nearly linear increase of intensity with length 
of light-path. The observed doublet ratio testifies that, if the num- 
ber of atoms is doubled by matching each atom with a similar addi- 
tional atom, the line intensity is increased by 20 per cent; on the 
other hand, if the number of atoms is doubled by extending the 
light-path twice as far into space, the line intensity is nearly doubled. 
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The properties of sodium gas cannot therefore be identical at all 
points in interstellar space. Motions corresponding to the standard 
concept of galactic rotation do not suffice to solve the enigma. The 
only reasonable solution which has occurred to us is that of discrete 
clouds of sodium vapor, although on account of the complexity of 
the problem we are scarcely disposed to argue that other possibilities 
are entirely precluded. Such discrete clouds perhaps form the gase- 
ous counterpart to the well-known dust clouds of obscuring matter 
along the Milky Way. The origin and stability of these clouds may 
furnish interesting questions for discussion. 
CARNEGIE INSTITUTION OF WASHINGTON 


Mount WILSON OBSERVATORY 
May 10937 
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PHYSICAL PROCESSES IN GASEOUS NEBULAE 
Il. THEORY OF THE BALMER DECREMENT 
DONALD H. MENZEL AND JAMES G. BAKER 


ABSTRACT 

An exact algebraic solution of the equations given in the first paper of this series is 
carried out for the statistical equilibrium of an assembly of hydrogenic atoms. The ex- 
citing radiation is assumed to come only from beyond the limit of the Lyman series. 
Hence discrete states are filled only by capture from the continuum and by cascade from 
higher discrete states. The electrons are assumed to have a Maxwellian velocity dis- 
tribution. The formula is given in terms of a quantity bn, the ratio between the actual 
population of state ” and the population under conditions of thermodynamic equilibri- 
um. An asymptotic expression is given for bn in the case of large quantum number, 
from which it is shown that bn approaches unity as m becomes infinite. Numerical solu- 
tions are deferred until the next paper. 

The problem of the theoretical decrement of the Balmer series has 
been attacked by several investigators, notably Plaskett,’ Carroll,’ 
and Cillié.* General equations for the statistical equilibrium of an 
assembly not in thermodynamic equilibrium were first given by 
Rosseland.‘ For each atomic state one writes down the condition 
that the number of transitions from that state, by all possible proc- 
esses, is equal to the number of transitions to that state, per unit 
time. Equilibrium is thus insured. Atoms possess, however, an in- 
finite number of states, and an accurate solution of the problem re- 
quires the treatment of an infinite number of simultaneous linear 
equations. 

To simplify the problem, Plaskett and Carroll have assumed that 
the atom possesses but seven discrete levels, the eighth denoting 
ionization. The effect of the remaining levels was entirely neglected. 
For the atomic emission probabilities they used the values of the 
Einstein A’s as given by Slack’ and of the capture probabilities as 
given by Gaunt.® Cillié?7 has recently employed the more extensive 
data on the emission probabilities, as given by Menzel and Pekeris,* 


to extend his original calculations to fourteen states. 


t Harvard Circ., No. 335, 1928. 5 Phys. Rev., 31, 527, 1928. 

2 M.N., 90, 588, 1930. 6 Phil. Trans., A, 229, 163, 1930. 
3 [bid., 92, 820, 1932. 7M.N., 96, 771, 1936. 

4Ap. J., 63, 218, 1926. 8 M.N., 96, 77, 1935. 
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One further simplifying assumption, justified to some extent by 
physical reasoning, is made. In the problem of the planetary nebula 
the argument is that the Lyman line radiation from the nuclear star 
is so depleted by the star’s atmosphere and the lower nebular layers 
that the only radiation available for atomic excitation lies beyond 
the limit of the Lyman series. 

For the present we shall accept the conditions of the preceding 
paragraph. It is instructive to carry out the accurate solution of the 
equations in this particular case. The general procedure followed 
by the earlier investigators is, however, open to question. The nu- 
merical solutions are none too satisfactory, because of the approxi- 
mations involved. One has the feeling that more attention should 
be paid to the neglected transitions. Continuity of the energy levels 
through the higher terms and on into the continuum is an essential 
feature of every atom. Furthermore, the use of the exact values of 
the Einstein A’s is incompatible with the use of the Kramer-Gaunt 
approximation to the emission in the continuum, since such pro- 
cedure involves a discontinuity at the series limit. A further and 
serious objection is that this method gives no information about the 
intensities of high series members. 

The exact equations for the statistical equilibrium have been given 
in the first paper? of this series, hereinafter referred to as “I.”’ An 
atom in quantum state m may make spontaneous transitions to low- 
er discrete states m’, or, under the influence of radiation, to either 
higher discrete states m’’ or higher continuous states x. Let Fz, de- 
note the number of transitions n to ’ per unit volume per second. 
Then the condition of equilibrium (I, equation 31) requires that 


io) foo] 
+ Fendv + = > Fan” 


n’=1 n”=n+1 
(1) 


+ + 


n'=1 


Under the assumed conditions, there is no radiational excitation 
from any quantum state but the first; and, since the effect of Lyman 


9 Menzel, A p. J., 85, 330, 1937. 
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line radiation is neglected, the only excitation of the 1-quantum 
atoms is to the continuum. Hence, for each state other than the 
first we have an equation of the form 


oo on n—1 
Fen f = Fan’. (2) 


n'=1 


Substituting from I, equations (10) and (25), we find that 


< 
where 
Sn = e*n—E(—X,)], 
mae hRZ? _ 1.5703 X 10% 
for hydrogen. 


In the foregoing equations NV; and JN, are, respectively, the num- 
ber of ions and electrons per unit volume; 7, is the electron tempera- 
ture; g, the Kramer-Gaunt factor for the transition in question; and 
b,, a factor that denotes the departure of the assembly from thermo- 
dynamic equilibrium at temperature 7.. It is defined by the rela- 
tion (I, 4) 


= 
Nn b,N iN. (ommkT.)3/2” ° (4) 
For thermodynamic equilibrium, 5, =1. 
— (5) 


K = mo 
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To duplicate Cillié’s solution we have merely to substitute the ac- 
curate g values for the first fourteen equations, cut off all the sum- 
mations at n”’=14, set Z (the averaged g for the bound-free transi- 
tions) equal to unity, and, neglecting all higher equations, solve 
simultaneously for the various 0,’s. We prefer to start from the op- 
posite end and obtain the b,’s for the higher levels first. 

A consistent procedure is to set, as a first approximation, g equal 
to unity throughout. The second of the two summations is easily 
evaluated for large values of n: 


n—I 


n? 2 I I 
| (6) 
I I 

[sin 39 In 2 (7) 
~3lnn+ 1.04. (8) 


For as small as 10 the bracket in (7) is correct to five significant 
figures. In deriving these formulae we have noted the identity 


(9) 


and also the well-known expression for the sum to a terms of the 
harmonic series 


a 


I 
(10) 


where ¥ is Euler’s constant, 0.577..... Various functional expan- 
sions of an obvious nature have also been used. 

For the middle term of equation (3), certain relationships may be 
derived. 
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For all values of X, different from zero, 


Sn = — E(—X,)] = 
Ietk: 
x3 (11) 
3° 3: 


For X,>>1, it is convenient to employ the following asymptotic 
expansion: 


Equation (12) is needed only for the lowest quantum states. 

Thus, the only troublesome term in equation (3) is the first. The 
difficulty arises from the occurrence of the factor b, under the 
summation sign. Before proceeding to the exact solution of the equa- 
tion, we shall employ an approximation, which indicates clearly the 
physical nature of b,. From the character of the summation we note 
that the principal contribution to its value comes from n”’ terms bor- 
dering on . If, therefore, 6,” is a slowly changing function of »’’, 
we may remove it from the summation sign in order to obtain the 
first approximation. 

The first term may be written as follows: 


n? 


since, when »”’ is large, we may approximate the sum by an integral. 
If we change variables by 


ARZ (1 I 


2hRZ? dn” (14) 


n''3 = dy 
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the summation becomes, by equation (11), 


—y 


(15) 


When the limiting formulae (8), (11), and (15) are substituted 
back into equation (3) and the common factors are canceled, we 
may solve for b, and find 


As n>, b,>1. This limiting value of unity is in accord with our 
physical intuition, which leads us to expect a continuity between 
states of negative and positive energy near the ionization limit of the 
atom.’® Since our calculations are based on the assumption that the 
electrons possess a Maxwellian velocity distribution, the value of }, 
appropriate to the continuum where thermodynamic equilibrium ob- 
tains, is unity. The approximate value of (16) can be employed, 
however, only for the uppermost states. 

Equation (3) may be written as an infinite set of simultaneous 
equations involving an infinite number of unknowns. One may write 
down symbolically the infinite determinants for the general case and 
expand by minors. Since equation (3) does not contain n’ directly, 
the solution for any b, is of precisely the same form. Accordingly, 
the equations as written for the solution of b, begin with the mth 
equation and not with m= 2, although the answer is the same. This 
fact simplifies the solution considerably. The determinant of the 


coefficients of the unknowns has the value t;. Since ¢;>1, this 


expression is, of course, infinite in value. The determinants in the 
numerator, however, contain the same minors from some point on, 
and hence, the ratios reduce to infinite series. 

10 The limiting form of bn was given by Menzel and Cillié, Ap. J., 85, 88, 1937. 
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An expansion of the terms and a recombination of the cyclical 
subscripts reduce to the exact solution 


Sn t Va 
b, = eXnt, (17) 
where 
Sittin 
n+I n+1 j+1 
> (18) 
4 3 > 
«tr 
It is evident that these summations converge individually, for 
and 
= 
n+1 
Indeed, since 
b,> 1, V,>lnn + const. 


An examination of equation (3) shows that formula (17) may be 
derived directly by the method of successive approximations. This 
procedure is by far the shortest. The fact that these two entirely 
different methods yield precisely the same formula for 6, is interest- 
ing in itself. What is so obvious here, however, will not be so at all 
for the more complicated solutions to come. Accordingly, we have 
employed the determinant method to indicate how the future re- 
finements will go. 

We are now enabled to determine precisely the numerical values of 
b,. The solution (17) is of the final form for the nebular problem, 
since, upon introducing the Gaunt factor and collision and radiation 
terms from I, we need only to redefine the functions given above. 
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The next paper will contain useful tables of the important func- 
tions of the preceding equations and, in particular, will contain a 
tabulation of the values of b, over a wide range of temperature and 
of quantum number, on the basis of the straight capture spectrum of 
hydrogenic atoms. A comparison with observations and results of 
other investigators will also be included. 


HARVARD COLLEGE OBSERVATORY 
May 14, 1937 
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THE OPACITY IN THE INTERIOR OF A STAR 
S. CHANDRASEKHAR 


ABSTRACT 
In this paper two integral theorems on the radiative equilibrium of a gaseous star 
are proved. 


I 

In two recent papers’ the author has proved some general theo- 
rems on the equilibrium of a star. These theorems are of some im- 
portance in the theory of stellar structures, in so far as they provide 
inequalities for the physical variables, e.g., central pressure, mean 
pressure, central radiation pressure (1 — .), etc., which should be 
valid under very general circumstances. The method consists in 
obtaining inequalities for the physical variables which are direct 
consequences of the equation of hydrostatic equilibrium: 


(1) 


dr r 


In obtaining inequalities based on equation (1), one generally re- 
stricts one’s self to such equilibrium configurations as are character- 
ized by the mean density f(r), inside r decreasing outward. Further, 
in obtaining inequalities for the mean temperature a further restric- 
tion, namely, that (1 — 8) decreases outward, is introduced (cf. II, 
Theorem 7). 

In this paper we shall obtain certain inequalities for equilibrium 
configurations in radiative equilibrium. We shall then have an addi- 
tional differential equation for the radiation pressure p,(= 3a7"), 
namely, 

dp, _ xL(r) (2) 


dr 4mcr? 


where x is the opacity coefficient, c is the velocity of light, and L(r) 
is the amount of energy crossing the spherical surface of radius r. 
The numbering of the theorems is continued from II, and refer- 
ences to the equations of that paper are inclosed in square brackets. 
*M.N., 96, 644, 1935; Ap. J., 85, 372, 1937. These papers will be referred to as 
“T” and “II,” respectively. 
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Il 


In the following, ZL, M, and R refer to the luminosity, the mass, 
and the radius of a star; and the auxiliary variable 7 is defined by 


(3) 


THEOREM 8.—In a wholly gaseous configuration in radiative equi- 
librium, in which the mean density p(r) inside r decreases outward, we 
have 

cGM(1 — 
(4) 


Kn 


where (1 — B*) has the same meaning as in Theorem 2 and kn is 
defined by 


R 
Part of the analysis leading up to this theorem is originally due to 


B. Strémgren.? 
Proof: The equation of radiative equilibrium (2) can be written as 


— 
dp, r4 (6) 
Since 
— & M(r)dM(r) 
dP r4 (7) 
we have, using equation (3), that 
L 
dp, = 4mcGM xndP (8) 


Integrating equation (8) and using the boundary condition that 
pr = oatr = R, we have 


L 
pr = xndP . (9) 


2 Handbuch der Astrophysik, 8, 159, 1936. 
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The foregoing equation has been given before by Strémgren. From 
equation (g) we can easily show that if xn decreases outward, (1 — B) 
would also decrease outward. Extending the integral from o to R and 
using the definition (5) for the average value for xn, we clearly have 
that 


4ncGM Kn, (10) 
or 
4mcGM (i — 


Since we have assumed that the mean density decreases outward, 
we can apply Theorem 2, which states that 

B*, (12) 

where #* satisfies a certain quartic equation (equation [3]) and is 


determined by the mass M, uniquely. Combining equations (11) and 
(12), we have 


* 
n 
which proves the theorem. 
Ill 


THEOREM 9.—In a wholly gaseous configuration in which the mean 
density p(r) inside r and the rate of generation of energy ¢€ decrease 
outward, we have 


(14) 


where x is the mean opacity coefficient defined by 


° 
pam xdP , (15) 
R 


and the equality sign in equation (14) is possible only when € is con- 
stant. 
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Proof: This is an immediate consequence of Theorem 8. For, if 
e decreases outward, 7 must also decrease outward, and consequent- 
ly the minimum value of 7 is unity. Hence 


kn 2 K, (16) 
the equality sign in equation (16) being possible only when n= 
constant = 1, i.e., when ¢ is constant. By Theorem 8 


kn < (17) 


Combining equations (16) and (17), we have the required result. 


IV 


We will apply equation (14) to certain practical cases of interest. 
Numerically, equation (14) reduces to 


K < 1.318 X 10! (1 — 6*), (18) 


where Lo refers to the luminosity of the sun. 

For Capella we have M = 4.18© and L = 126 Lo. Assuming 
bh = 1, the solution of the quartic equation for B* yields 1 — B* = 
0.22. Hence 


KCapella < 96-1 gm" cm? . (19) 
In the same way for the sun, we find (u = 1, 1 — B* = 0.03) 
Ko < 395 gm" cm’. (20) 


V 
There is one interesting application of equation (14) to stellar 
models in which the opacity coefficient x is assumed to be constant. 
For equation (14) can then be written as 
4mcGM (1 — B*) 


K 


(21) 


L* 


Equation (21) has to be interpreted in the following sense: If 
L > L*, then the configuration must be characterized by the oc- 
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currence of negative density gradients (i.e., p(r) increases outward 
in some finite regions of the interior), no matter what the law of 
energy generation is, provided only the rate of generation of energy 
e decreases outward. On the other hand, if L < L*, it does not 
necessarily follow that the configuration is characterized by positive 
density gradients throughout its interior. But if L < L* we can al- 
ways find a “‘mild’’-enough law for the rate of generation of energy 
such that the configuration is characterized by a positive density 
gradient throughout its interior. Further, it should be noticed that 
if 


L> isn (22) 


then no equilibrium configuration is possible. The inequality (22) 
is interpreted by the statement that if Z > L, then the configuration 
would “blow up.” We now see that this tendency to “blow up” 
must set in at lower values for the luminosity, in the event of nega- 
tive density gradients in its interior. Negative density gradients 
must certainly exist for configurations with L > L*. Depending on 
the concentration of the energy sources toward the center, the nega- 
tive density gradients will set in for some L < L*. 


VI 


Corollary to Theorem 9.—If, in addition to the conditions of Theo- 
rem 9, x is assumed to increase outward, then 


(23) 


where x, is the opacity at the center. This is, of course, obvious. 

If we assume any definite law for opacity, then equation (23) can 
be converted into an inequality for the central temperature, for a 
star of known mass and luminosity. Thus, if we assumed that 


(S>0), (24) 


3E. A. Milne, M.N., 91, 4, 1930. See esp. pp. 12, 13, and 53 of this paper. 
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then x would increase outward if (1 — 8), and T would decrease 
outward; for we can write equation (24) as 


(25) 


(26) 


Combining this with (23), we have 


Lk, a 
4ncGM k 3 (1 — B*)?° (27) 
The foregoing inequality giving the minimum central temperature 
for a star of known M and Lis generally not as good as the minimum 
central temperature set by Theorem 7 for a star of known M and R. 


APPENDIX 


We have seen that (1 — 8*), giving the maximum possible (1 — 6.) in a 
wholly gaseous configuration in which the mean density f(r), inside r, is as- 
sumed to decrease outward, plays an important role in Theorems 7, 8, and 9. 
It is therefore convenient to have a table giving M for different values of 
(1 — 6*). Table I should be sufficient for most purposes. 


TABLE 1 
SOLUTIONS OF THE EQUATION 


YERKES OBSERVATORY 
June 8, 1937 
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THE ANNUAL VARIATION IN THE INTENSITY 
OF THE ZODIACAL LI ;HT* 


C. T. ELVEY 


ABSTRACT 


The series of visual observations of the zodiacal light made by the Japanese observers 
has been reduced to a uniform system of intensities with the aid of photoelectric ob- 
servations of the standard comparison areas in the Milky Way. Corrections were ap- 
plied for differences in the zenith distance of the zodiacal light and the comparison field. 

Annual variations are shown to exist for both the evening and the morning zodiacal 
lights. The maximum for the evening light is near the middle of January, and that for 
the morning light is near the end of the preceding October. 

The suggestion is made that the enhancement of the zodiacal light at certain seasons 
of Ac year may be due to scattering by meteoric material associated with cometary 
orbits. 


In a recent survey of the light from the night sky with a recording 
photoelectric photometer, Elvey and Roach’ obtained the intensities 
of the zodiacal light over the sky. The observations were grouped 
roughly in accordance with the seasons, and it was shown that there 
is a seasonal variation in the intensity of the zodiacal light during 
the period covered by the observations, namely, November, 1934, to 
July, 1935. The maximum of intensity of the morning and evening 
zodiacal light occurred, in general, at the times of the best visibility 
of these phenomena. To account for the displacement of the maxi- 
mum of the morning light with respect to that of the evening light 
it was suggested that there may be a condensation in the scattering 
medium, which is in conjunction with the sun in the late winter. 

Since the above-mentioned photometric observations are some- 
what limited in number and time, they are probably affected by 
variations in the intensity of the zodiacal light which are irregular in 
character and which may or may not be due to differences in the 
transparency of the atmosphere. For this reason it seems advisable 
to investigate the large mass of visual observations for possible 
seasonal variations, even though they do not have the accuracy of 
photometric data. An excellent series of observations were obtained 
by the Zodiacal Light Section of the Society of Astronomical Friends 

* Contributions from the McDonald Observatory, University of Texas, No. 4. 

J., 85, 213, 1937- 


84 


THE ZODIACAL LIGHT 85 


of Kyoto, Japan, which later became the Oriental Astronomical 
Society. These observers, numbering as many as twenty or more, 
have been making systematic observations of the position and out- 
line of the zodiacal light as well as of its color and brightness and of 
such other miscellaneous observations as pulsations in brightness 
and variations in outline. 

We are concerned only with the estimates of the brightness of the 
zodiacal light, which we have compared with selected regions of the 
Milky Way. The Japanese reports do not give details of the posi- 
tions of the zodiacal light observed and that of the comparison field; 
however, approximate values can be obtained from the statement 
that the maximum of the zodiacal light was observed. Japanese 
Central Standard Time is given for each observation. From these 
data we see that the observations have been taken over a consider- 
able range in relative zenith distance between the zodiacal light and 
the comparison field, which in many cases gives systematic effects. 
For example, a given comparison field, such as that in Lacerta and 
Cepheus, was used over a period of two or three months in succes- 
sion. The zodiacal light remains at nearly the same zenith distance 
throughout the series, but the comparison field changes month by 
month, thus introducing into the observations an effect which may 
be seasonal. 

Before discussing further the calibration of the observations or 
the corrections to be applied, we should look, perhaps, at the quan- 
tities which are measured by the eye. The eye measures the intensity 
of an object relative to the intensity of the background. In the case 
of the comparison fields this is the integrated light of the faint stars 
of the Milky Way, plus the scattered light, plus the galactic light, 
while the background illumination consists of the faint stars, the 
auroral light, and the general amount of zodiacal light over that re- 
gion, as well as the scattered light from all these sources. The light 
which originates outside the atmosphere will be affected by extinc- 
tion, which is a function of the secant of the zenith distance up to 
values of about 70°, when some other function of the zenith distance 
must be used, such as that derived by Bemporad.? The light origi- 
nating in the earth’s atmosphere increases in intensity with zenith 


2 Handbuch d. Astrophysik, 2, 184, 1929. 
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distance, except for very large values, when the extinction becomes 
powerful enough to cause a reversal and, consequently, a maximum 
of intensity at about 80° zenith distance. Since the comparison 
field, located in the Milky Way, is outside the atmosphere and con- 
sequently decreases in intensity with zenith distance by extinction, 
while the background is composed mainly of light originating in the 
atmosphere which increases with zenith distance, it is readily seen 
that extinction and variation of the background combine to cause 
a decrease in the contrast between the Milky Way field and the 
background. 

In the case of the zodiacal light the observations are made at the 
point of maximum intensity. This point, except in rather adverse 


TABLE 1 
Intensity in Numbers of Stars 
Comparison Field of the Tenth Magnitude 
per Square Degree 
300 
200 


weather conditions, should be near the same zenith distance, on the 
average. The position is dependent upon the extinction and back- 
ground illumination; and since we do not have data as to the altitude 
at which the zodiacal light was observed, I have estimated it as 20° 
(zenith distance 70°). 

In order to convert the estimates of the intensity of the zodiacal 
light into true intensities, it is necessary to know the intensities of 
the comparison fields. These should be determined with a visual 
photometer; but since such measures are not available to me, I have 
used measures made with a photoelectric photometer. For a study 
of a variation in the zodiacal light the use of photoelectric intensities 
is satisfactory, unless there are differences in color between the vari- 
ous comparison fields. The intensities of the comparison fields are 
the averages over the Milky Way in the designated positions and 
are listed in Table 1 to the nearest 25 units. 


| 


THE ZODIACAL LIGHT 87 


All estimates of the intensity of the zodiacal light for the years 
1928-1935, inclusive, published in the Kyoto Bulletins, were re- 
duced to intensities by means of the foregoing calibrations. It was 
soon evident that there were curious discrepancies within the ob- 
servations over a short period of time—say during one month. 
These are related to the comparison field employed by the observer. 
In cases of large differences it was found that one comparison field 
was usually at the same zenith distance as the zodiacal light, while 
the other was at a small zenith distance, thus indicating that extinc- 
tion and background illumination had an important effect on the ob- 
servations. 

It is necessary to know the amount of the sky background at each 
point of observation and also the value of the extinction, in order to 
reduce the observations to a uniform system of intensities. In the 
absence of such data for the place of observation, we can only try 
to estimate a reasonable value for the two quantities from our experi- 
ences at other places. The series of photoelectric observations made 
by Elvey and Roach’ was used for this purpose, and it was found that 
the average extinction was 0.25 mag. at the zenith. This is for the 
blue end of the spectrum. Of course, the extinction would be smaller 
in the visual region. On the other hand, our observations were made 
at an elevation of about 2000 meters, and the Japanese observations 
of the zodiacal light were made close to sea-level, hence the extinc- 
tion would have to be increased. The foregoing value of 0.25 mag. 
has been assumed for the extinction to be used in correcting the 
zodiacal light observations. This is quite close to the extinction 
constant, 0.23 mag., used in the mean extinction tables for visual 
light at Potsdam,‘ at 100 meters elevation, and to the value of 0.195 
mag. in Bemporad’s? extinction table for sea-level. A table of differ- 
ential extinction corrections can be readily computed, either in mag- 
nitudes or as a factor, to correct the observed intensities. The differ- 
ential extinction corrections are listed in columns 2 and 3 of Table 2. 

A correction for the variation in the background intensity with 

3 Bull. Kwasan Obs., Nos. 142, 164, 166, 168, 179, 190, 208, 212, 252, 255, 273, 275, 


299, 325 (1928-36). 
4 Handbuch d. Astrophysik, 2, 264, 1929. 
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zenith distance was taken also from the observations with the photo- 
electric photometer. The average intensity of the auroral light plus 


TABLE 2 
CORRECTION FOR AURORAL 
CoRRECTION FOR Extinction | LicHt PLus SCATTERED LIGHT T 
ZENITH FOR AREA IN AURIGA pai 
CoRRECTION 
DISTANCE 

(FACTOR) 

Mag. Factor Mag. Factor 
+0.47 1.54 +o0.22 1.22 1.89 
eae + .47 1.54 + .20 1.20 1.85 
+ .45 1.51 + .17 1.77 
+ .44 1.50 + .16 1.16 
+ .43 1.49 + .15 1.71 
+ .42 1.47 + .13 1.13 1.66 
+ .40 1.45 + .II 1.60 
+ .37 at + .08 1.08 
+ .33 1.36 + .07 1.06 45 
+ .28 1.30 + .05 1.05 1.36 
+ .22 1.22 + .03 1.03 1.26 
+ .13 + .02 1.02 
.00 1.00 .00 1.00 1.00 
— a4 0.80 .03 0.97 0.78 
—0.68 0.53 +0.03 1.03 0.55 

220 T T T T T T T T 
200+ + 


70 80 90° 


Fic. 1.—The average variation of auroral light plus scattered light, with zenith 
distance. 


scattered light, as observed at the McDonald Observatory, was used, 
and its variation with zenith distance is shown in Figure 1. The in- 
tensities are expressed in numbers of tenth-magnitude stars per 
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square degree. Since this observed curve would have to be corrected 
to apply to visual light at sea-level, and since we do not know the 
amount of the correction, these data have been assumed. Some of 
the changes for the variation in elevation would result in an increase 
of the intensities, while others would decrease them, thus justifying 
the use of the foregoing data as they stand. 

By assuming the curve in Figure 1 as representing the variation of 
the sky background with zenith distance, we can compute the ap- 
parent variation in intensity of a constant amount of illumination 
superimposed on the sky background. For purposes of illustration 
an intensity of 1oo has been added to the background, and the 
amount by which it is brighter than the background has been com- 
puted in magnitudes. This was converted into the factor for cor- 
recting the observed intensities, and these data are tabulated in 
columns 4 and 5 of Table 2. The total correcting factor is given in 
column 6, this being the correction for the comparison field in Auriga. 

It is seen that these correcting factors are large. It is, therefore, 
important that observations be made of points in the zodiacal light 
whose zenith distances are the same as those of the comparison field 
in the Milky Way. If the visual observers will take such precautions, 
this will add very materially to the value of their estimates of the in- 
tensity of the zodiacal light. 

The observations of the zodiacal light were not corrected individ- 
ually for the differences in zenith distance but were first averaged 
for a given month. The average zenith distance of the comparison 
field’ for the month was readily obtained, and the correction was 
then applied to the mean intensity for the month. The accuracy of 
the method and of the visual estimates hardly warrants the greater 
amount of work required in reducing each observation separately. 

The results are shown in Figure 2 in the form of curves. There 
were no obvious periodic variations in the observations, other than 
the annual variation; hence the data were arranged into annual 
curves for the evening and the morning zodiacal light. The total 
number of observations for each month is indicated by the size of the 
plotted point. The observations are plotted for the middle of the 


5 Usually only one comparison field was used over a period of one month. Otherwise 
it is necessary to make corrections for each comparison field of a different intensity. 
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month. A continuous heavy curve has been drawn for all the means 
containing twenty-five or more observations. During the summer, 
when the ecliptic is low in the southern sky, there are very few ob- 
servations and the means have low weight. 


4 
700 
600+ 4 
500} 
oo} 
4 \ ra 
\ ¢ \ 
100} 
Jan. 15 July 15 Jan. 15 


Fic. 2.—The annual variation of the zodiacal light. The upper curve is for the evening 
observations, and the lower for the morning observations. The size of each dot repre- 
sents the number of observations forming the mean, the large dots being for more 
than one hundred observations and the small circles for less than five. 


In both the evening and the morning zodiacal light there is a 
rather sharp maximum, the one for the evening light falling very 
near the middle of January and the one for the morning light near 
the end of the preceding October. There is a broad minimum extend- 
ing through the season when it is difficult to observe the zodiacal 
light, namely, from May to October for the evening light, and from 
January to June for the morning light. 
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The time of maximum intensity of the morning zodiacal light 
agrees satisfactorily with that found by Elvey and Roach from their 
photoelectric observations. The maximum intensity of the evening 
zodiacal light does not agree: the photoelectric observations placed 
it in April, instead of in January. It is probable that because of the 
small number of photoelectric observations the true seasonal varia- 
tion was masked by rather large irregular variations in the intensity 
of the zodiacal light. 

No mention has been made of a correction for a variation in the 
elongation of the observed point in the zodiacal light through the 
year. When the observations are made with the ecliptic in its most 
favorable position, the point is much closer to the sun than in its less 
favorable positions. We can obtain an estimate of the difference in 
longitude of the point observed and the sun by taking the average 
of the times of the observations and the assumed zenith distance of 
the observed point. For the observations in the more favorable posi- 
tion of the ecliptic the differential longitude is between 40° and 45°. 
For the summer observations, or the times when the ecliptic is in- 
clined toward the horizon, the differential longitude is estimated at 
about 70°. By referring to the isophotes of the zodiacal light deter- 
mined by Elvey and Roach,’ it is estimated that a factor of 2 or 2.5 
must be applied to the observations of the zodiacal light made during 
the unfavorable season in order to reduce them to the same units as 
the observations taken when the ecliptic is nearly perpendicular to 
the horizon. In Figure 2 the probable intensity of the zodiacal light 
is indicated by the faint continuous line. 

Attention should be called to one discordant point in the morning 
zodiacal light. The observations for July indicate a large value for 
the intensity of the zodiacal light during that month. The plotted 
point is an average of the observations for only two years; in addi- 
tion, all of the observations for one year were taken under bad ob- 
serving conditions. The observations for the other year, however, 
are apparently good. It will be rather difficult to decide from ob- 
servations taken only at latitudes greater than 30° north whether the 
foregoing discordant point is real. We do not have available photo- 
metric observations of the zodiacal light taken from the southern 
hemisphere, but there are some sketches of the zodiacal cones made 
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by R. B. Bousfield® at Stanthorpe, Australia, which indicate that 
during July the morning zodiacal light reaches large intensities. 
The observations of the intensity of the zodiacal light, after due 
consideration has been given to the accuracy of the observations 
and to the corrections which had to be applied to reduce the inten- 
sities to a homogeneous system, indicate rather definitely that the 


180° 


Fic. 3.—Positions of the earth in its orbit at the times of maximum intensity of the 
zodiacal light, with the directions toward the points of observation, in elongations of 
40°. The enhancement of the zodiacal light may be a result of meteoric material moving 
in a cometary orbit. Three possible orbits are shown: A, Daniel’s Comet; B, Comet 
1813 I; and C, Comet 1819 IV. , 


annual variation noted is real. However, no interpretation of the 
annual variation of the zodiacal light can be definitely made without 
an investigation of the distribution of its intensity over the celestial 
sphere in conjunction with the theories of scattering of light by 
particles. Besides the photometric observations giving the distribu- 
tion of intensities, a study of the color of the zodiacal light will give 
additional information concerning the size of the particles producing 
the scattering of sunlight. 


6 M.N., 94, 824, 1934. 
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For the present we can say that the observed maxima represent a 
greater density of the scattering particles in the directions of the 
observations, namely, at an elongation of 40° west of the sun for 
the morning light and 40° east of the sun for the evening light. This 
is represented graphically in Figure 3, which shows the positions of 
the earth on the two dates of maximum intensity of the zodiacal 
light and the directions of the observed points. Data are not avail- 
able for other points, and a limit for the cloud of particles cannot be 
given. However, the fact that the maxima are rather sharp indicates 
that the cloud of particles is not very extensive. 

The maxima are present each year and at the same time; it is, of 
course, very improbable that there is a stationary cloud of particles 
at the indicated position, in longitude 77° and at 0.65 astronomical 
units from the sun. The observed effect could perhaps be produced 
by a meteor swarm seen lengthwise along a considerable part of its 
orbit. However, none of the observed orbits of meteor swarms fits 
our data. There remains the remote possibility that one of the known 
comets has produced a meteor swarm in the required position. An 
examination of Crommelin’s Comet Catalog shows three comets that 
would come under consideration. The orbits of these comets are 
shown in Figure 3. Daniel’s comet, 1907 IV, is probably too close to 
the sun. The next orbit belongs to Comet 1813 I, and the outer one 
is for Comet 1819 IV. It is, of course, not possible to prove this con- 
nection with the zodiacal light, and the association with cometary 
orbits can only be tentatively suggested. Those indicated are either 
parabolic or have very long periods, and it is doubtful whether such 
comets could shed sufficient meteoric material to form a continuous 
swarm of particles around the orbits. 
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| AN EMISSION NEBULA NEAR o SCORPII* 
OTTO STRUVE 


ABSTRACT 
With the use of a photographic emulsion sensitized for Ha, diffuse galactic nebulae 
have been classified into three groups, depending upon the relative intensity of the red 
image to the violet image. This method is used to determine the character of nebular 
spectra. A nebula near o Scorpii is found to be of the emission type. The red image 
does not, however, coincide exactly with the violet image. A part of the latter is there- 
fore attributed to pure reflection. The spectrum of the star is Br. | 
Photographs of the nebulous region in Scorpius and Ophiuchus, | 
obtained last year with the 94-mm Schmidt camera, revealed’ that 
the bright nebulae surrounding the B-type stars 22 Scorpii (B3), 
p Ophiuchi —24°12684 (B3), o Scorpii (Br), and Scorpii 
(B3) are strong in violet light. This was interpreted to mean that 
the nebulae were at least as blue as, and possibly slightly bluer than, 
.the illuminating stars. However, this interpretation stands or falls 
; with the assumption that the spectra of the nebulae are continuous. 
WW The spectra of the nebulae near v Scorpii, 22 Scorpii, and p Ophiu- 
chi have actually been observed as continuous.? The nebula near 
a Scorpii, however, has not been observed, and Hubble has merely 
inferred from the spectral type of the star, B1, that the spectrum of | 
the nebula is probably continuous. 

This nebula is so faint, and the star near it so bright, that it is 
exceedingly difficult to obtain a spectrogram which would definitely 
settle the question. However, following a suggestion by Keenan, 
Elvey and I successfully used red-sensitive plates with red filters to 
isolate the Ha-radiation.4 Emission nebulae are very strong when 
photographed on such plates. 

This method has now been materially improved by the use of a 
new emulsion furnished by the Eastman Kodak Research Labora- 


* Contributions from the McDonald Observatory, University of Texas, No. 5. 

t Ap. J., 84, 219, 1936. 

2 Hubble, ibid., 56, 405, 1922. 

3 Ibid., 84, 604, 1936. 4 Pub. A.A.S., 9, 23, 1936; Ap. J., 85, 252, 1937. 
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tory. This emulsion‘ is especially sensitive to the spectral region on 
both sides of Ha. It is much slower than emulsion IC for yellow 
light but surpasses it at Ha. The contrast is very good. 

Plate I shows the region near y Cygni. Exposure (a) was ob- 
tained on an Eastman Process emulsion, without filter. This repre- 
sents the nebulae in violet light, since the Process emulsion is sensi- 
tive only to about \ 4800. Exposure (b) was taken on emulsion IG 
with a yellow filter, and exposure (c) on the Ha emulsion with a red 
filter (Wratten 29 F). This combination cuts off sharply at \ 60c0 
on the violet side and just beyond Ha on the red. It therefore trans- 
mits effectively only a relatively small amount of continuous radia- 
tion from the night sky, while any existing nebular emission in Ha 
is transmitted freely. Accordingly, the contrast between an emission 
nebula and the sky is enhanced. 

The nebulae near y Cygni are well shown on the violet exposure. 
They are almost completely absent in yellow light, but they are ex- 
ceedingly strong in Ha. This confirms our previous result® that these 
nebulae are of the emission type. The emission which completely 
surrounds y Cygni (near the center of the plate), being strongest east 
of the star, shows rather conclusively that y Cygni, type F8p, is re- 
sponsible for the nebular emission. Also of great interest is the faint 
film of nebulosity which covers the entire western half of the plate 
and which is obviously also due to emission. This sheet of nebulosity 
covers a considerable portion of the eastern and northeastern part 
of the great Cygnus cloud—a region which is rich in B- and O-type 
stars.’ 

Plate II shows two photographs of the region near o Scorpii. 
a Scorpii with its red nebula is on the left. The Ha emulsion was 
used in both cases, and the filter was Wratten 29 F. The blue nebula 
near 22 Scorpii is completely absent, in spite of the great contrast 
of the emulsion. The slightly redder nebulosity near — 24°12684 is 
barely visible. 

Quite unexpected is the presence of a large unsymmetrical nebula 
near o Scorpii. Violet-sensitive plates show a large, fairly symmetri- 

5 Emulsion Spl. Exptl. H-alpha, No. 63147 (Eastman Kodak Company, Rochester, 
New York). 

6 Ap. J., 85, 252, 1937. 7A.N., 231, 17, 192. 


| | 
| 
| 


96 OTTO STRUVE 


cal nebula. Barnard’s photographs’ show it slightly stronger on the 
east than on the west, and considerably stronger on the north than 
on the south. Ross’s photographs and those taken with the Schmidt 
camera’ confirm this appearance. On the Ha plate the nebula is 
fairly strong. However, it is almost entirely concentrated in one 
quadrant—from the north to the west, with a slight extension to- 
ward southwest but with only the barest trace east of the star. 

It is safe to infer that the nebula shown on the Ha emulsion has 
an emission spectrum; it is also certain that only a part of the light 
observed on the violet-sensitive plates can be identified with the 
emission nebula. The rest must be of the reflection type. Why only 
the western portion of the nebula emits Ha cannot be answered. The 
densest part of the dark nebula illuminated by o Scorpii lies north- 
east of the star. There is an arm of obscuring nebulosity on the west. 
The star density is greatest south and southeast of o Scorpii. 

The spectral type of a Scorpii is B1. The lines are fairly sharp, 
and there is not much uncertainty in assigning the proper class. 
The hydrogen lines are rather weak, suggesting an early subdivision 
of class B, but the Sz 1v lines \ 4089 and A 4116 are not strong enough 
to make it earlier than B1 and not weak enough to make it Bz. The 
existence of a mixed spectrum in connection with a Br star is not 
very unusual: Hubble lists at least one other case.? The unusual 
thing about the emission nebula near @ Scorpii is the fact that it does 
not coincide with the violet reflection nebula. 

A number of other nebulae have been observed with red-sensitive 
plates in order to establish the character of their spectra. The re- 
sults for all nebulae are summarized in Table 1. The last column 
contains a classification of the nebular image in red light. The desig- 
nations I, II, and III are relative to the image obtained on a violet- 
sensitive plate of high contrast. Comparisons were made either with 
Process plates of the same region taken with the Schmidt camera or 
with the photographs in the atlases by Barnard and by Ross. Ex- 
perience has shown that on the small scale of the Schmidt camera a 
Process plate reveals approximately the same amount of detail as is 
shown on the photographs of either atlas. Class I indicates that 


§ Atlas, Pl. 13. 
9 Ap. J., 56, 184, 1922. Star Bond 734 in NGC 1982, type Bip. 
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TABLE 1 
NEBULAE OBSERVED IN RED LIGHT* 
Exciting Sp. of Sp. of 
No. Nebula as Star Nebula Class Remarks 
IC 59 Cas Bop E Il or Ill I 
IC 63 y Cas Bop E Il 2 
Pleiades Merope Bs C Ill 3 
NGC 1499 Per Oes E I 4 
IC 2118 Rigel B8p Ill 5 
IC 434 E I 7 
NGC 2023 — 2°1345 B2 Il 8 
NGC 2024 ¢ Ori Bo I 9 
NGC 2175 +20°1284 I 9 
NGC 2237 Cluster Oes5 I II 
NGC 2245 +10°1159 Bip Ill 12 
S Mon S Mon II 13 
Large nebula | S Mon Il 14 
Sco a Sco B3 C  |land 15 
See IC 4604 p Oph B3 C Ill 16 
— 24°12684 — 24°12684 | Ill 16 
IC 4605 22 Sco B3 C Ill 16 
a Sco a Sco I 17 
+40°4103 Cyg? =e I 18 
Large nebula | y Cyg ? I 19 
Ane tee NGC 7000 a Cyg A3p C+E I 20 
57 Cyg a Cyg Sere I 20 
Cyg Cyg | II 21 


* Class I: the red images are stronger than the violet images. Class II: the red images are equal to 
the violet images. Class III: the red images are weaker than the violet images. The spectra of the nebulae 
are from Hubble (A. J., 56, 1922). 


REMARKS 
1. North of y Cas. 
2. East of y Cas. 
3. The nebula near Maia appears to be slightly bluer than that near Merope. 
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4. This nebula is not as strong in Ha as would be expected from the violet image. 

5. Large nebula northwest of Rigel. 

6. Outer nebulosity of Orion, northeast of NGC 2068. 

7. Bay nebula, south of ¢ Orionis. Not quite as strong in red light as NGC 2024. 
The extreme southern extension of this nebula seems to be relatively weaker in red light 
than the more intense northern portion. 

8. This nebula is unexpectedly red for the spectral class of its star. 

g. Very strong in Ha. 

10. The Ha plate is underexposed; and since the violet image is weak, the color 


class may not be correct. 
11. The Ha exposure is very strong; it shows greater extent and more detail than 


Ross’s Atlas, No. 35, or Lick Pub., 11, Pl. 27. 

12. The designation is according to Barnard (Pub. Lick Obs., 11, Pl. 28, 1913). 
Hubble designates this nebula IC 447 (Ap. J., 56, 405, 1922) and reserves the desig- 
nation NGC 2245 for a cometary nebula near by (ibid., 44, 195, 1916). 

13. Fora nebula stimulated by an Oes star the Ha image is surprisingly weak. It is 
too strong, however, to be attributed wholly to reflection. 

14. This is the faint diffuse nebulosity extending a degree or more northeast and 


east of S Mon. 
15. The eastern part of this nebula is of class III; the northwestern part is of class I. 
16. IC 4605 is bluest, and the nebula around — 24°12684 is slightly less blue. IC 


4604 is intermediate. 
17. Red reflection nebula near Antares. The strength of the yellow image proves 


that the nebula is not of the emission type. 
18. The nebula contains the Bzp star +40°4103; but the uniform surface brightness 
of the nebula, exhibiting no concentration toward this star, suggests that the excitation 


does not originate in it. 

19. This is the faint nebulous sheet extending toward the southwest from y Cyg. 
It covers a part of the cluster of B stars in the Cygnus cloud, 

20. The nebula between 56 and 57 Cyg is relatively stronger in Ha than NGC 7000 


(the North America nebula). 
21. The red star & Cyg is surrounded by a faint wisp of nebulosity which is strongest 


northeast of the star. It is probably of the reflection type. 


the Ha image is appreciably stronger than the violet image; II 
stands for approximately equal images; III is used to designate 
nebulae which are much weaker in red light than in violet light. 
Some of the red images used for this investigation were taken on 
emulsion IC with a red filter. 

Since many of the nebulae in Table 1 have known spectra, we can 
roughly calibrate our classification. Emission nebulae are almost al- 
ways lor II. The difference between these two classes is appreciable. 
Thus, the nebula near 15 S Monocerotis is relatively much weaker 
in red light than NGC 2237. Pure reflection nebulae illuminated by 
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stars of early spectral type are always of class III. Those illuminated 
by very red stars may be of class I or class II. 

The method described here has two advantages over the spectro- 
graphic method: it can be used for very faint nebulae, and it per- 
mits us to study the entire nebula—not only its brightest portion. 
Its principal disadvantage lies in the fact that it does not permit us to 
determine for each nebula the fraction of light concentrated in emis- 
sion lines and in the continuous spectrum. For example, we do not 
know whether the difference between emission nebulae of classes I 
and II is caused by a difference in the relative intensities of the 
emission lines or by different amounts of continuous spectrum super- 
posed over the emission lines. 


I am indebted to Dr. C. T. Elvey and to Dr. P. Rudnick for valu- 
able assistance in making the observations. 


YERKES OBSERVATORY 
AND 
McDonaLp OBSERVATORY 
May 19, 1937 
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NOTE ON A POSSIBLE EFFECT OF THE RECENT 
INCREASE IN BRIGHTNESS OF y CASSIOPEIAE 


ABSTRACT 


Attention is called to possible effects of the recent increase in brightness of y Cas- 
siopeiae on the brightness of the two near-by diffuse nebulosities. 


The rise in brightness of the Be star y Cassiopeiae in the latter 
part of 1936 may prove to be of additional interest because of pos- 


Fic. 1.—A is the projected distance in light-years from y Cassiopeiae (S) to the 
nearer nebula, N. The true distance is approximately A sec 6; the lag in time, 
observed at the earth £, between the stellar outburst and its effect on the nebula is 
approximately A (sec @+tan @). 
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sible effects on the two diffuse nebulosities which are situated close 
to the star. There can be little doubt that the nebulae owe their 
luminosity to the radiation of y Cassiopeiae, and it seems likely that 
the approximate doubling in the light of the latter would result in an 
appreciable increase in the brightness of the nebulae themselves. 

From similarity with Be stars observed by Trumpler’ in the 
cluster x Persei and in the Pleiades we can assume that the photo- 
graphic absolute magnitude of y Cassiopeiae probably lies between 
—1mag. and —3 mag. On the assumption that its apparent photo- 
graphic magnitude is 2, its distance from the earth lies between 40 
and 100 parsecs. 

If there is no foreshortening, the light of the star will reach the 
nebula (Fig. 1) A years later and will also be observed at the earth 
A years after the original outburst. If foreshortening is present, the 
true direction of the nebula makes an angle 6 with A, and the light 
from the nebula will reach the earth approximately 


A(sec + tan 6) [+ 9° > 6 > — 90°] 


years later than the observed rise in y Cassiopeiae. 

From the apparent observed shapes of the two nebulae it seems 
rather unlikely that @ > 45°. As the angular distance from the ver- 
tex of the nearer nebula to the star is 20’, the elapsed times at the 
earth between the observed rise in the star’s brightness and a pos- 
sible effect in the nebula are: 


0 x= 40 Parsecs | ~=100 Parsecs 
0.7 1.9 


Accurate photometric measures of the nebulae during the next 
few years may therefore show an increase in their brightness as a 
result of the increased strength of the radiation of y Cassiopeiae. 

W. W. Morcan 


YERKES OBSERVATORY 
April 21, 1937 


' Pub. A.S.P., 38, 350, 1926. 
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NOTE ON THE DUPLICITY OF NOVA HERCULIS 
ABSTRACT 


Micrometer measures and the visual appearance of the components of the nova are 
discussed, together with some conclusions derived therefrom. 

After the writer had found the nova to be a close double star in 
the 36-inch telescope of the Lick Observatory, micrometer measures 
were published by him' and by Votte,? Van Biesbroeck,’ Aitken,‘ 
and Simonov.’ The last author states that his images were only 
slightly elongated, and his distances are accordingly not reliable.° 
His results have not been used. 

The normal points representing these observations are given in 
Table 1. The two recent averages by the writer are from unpub- 


TABLE 1 
1935.515.......130°96 o”212 Am=o0.4 4 nK_ 36-inch 
128.0 .168 .5 5 VB 40-inch 
......... 131.5 .224 .4 5 K_ 36-inch 
2006. 131.6 .904 — .2 m Vou 24-inch 
128.3 .194 .6 5 m VB 4o-inch 
532.2 .505 .4 3,2" 4o-inch 
136.9 .39 .5 «1 36-inch 
1937.902....... 131.7 .62 .I 3 mK  4o-inch 


lished measures, which are in the table below. The observations 
made with seeing 3 are given half-weight. The angles and distances 


1936.687.......131°96 Am=o0.3 40, 1340 X seeing 4 
40, 1340 X 3 
1936.799....... 130.2 0.51 ©.4 40, 960 X 4 
13I.I 0.505 0.4 3,2n 
135°2 0°70 Am=o0.4 40,700 X_ seeing 3 
6.9 .© 40, 960, 1340 X 4 
131.0 0.65 ©.0 40, 960 X 4 

t Pub. A.S.P., 47, 228 and 267, 1935. 4 Pub. A.S.P., 48, 340, 1936. 

2 Ibid., p. 270. 5 A.N., 259, 81, 1936. 


3 Ap. J., 82, 433, 1935. ° Ibid., 256, 300, 1935. 
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in Table 1 are represented in Figure 1. My measures make the an- 
gles sensibly constant. The difference of 3° between these and Van 
Biesbroeck’s is certainly small for a pair of only 072 separation. Dr. 
Aitken’s result is, unfortunately, only a single night’s measure, with 
“seeing not of the very best, .... at least 3.” 

The distances appear to increase linearly with the time. I have 
drawn in the diagram a straight line which, in my opinion, is the 


T T T T T T T T T 7 


ESS: fr 


1935.0 1936.0 1937-0 


Fic. 1 


best that may be derived from the measures. It intersects the time 
axis at about 1934.82, which precedes the outburst on December 127 
= 1934.94 by 0.12 years. But this small difference is not considered 
significant. If it is possible to follow the nova for at least two more 
years, the reality of this difference may be discussed to advantage. 
The explanation of a difference of this nature would obviously be the 
gravitational deceleration of the two masses. 

Many papers have discussed models which would explain the ob- 
served spectral features. For mathematical reasons a spherical model 
is attractive and has been frequently used. But the visual observa- 


7 For instance, A. Beer, M.N., 97, 242, 1936. 
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tions show this to be very far from the truth. We shall discuss the 
visual observations of the two components, without trying to inter- 
pret them, since it is found that more data are needed. 

On July 3, 1935, when the duplicity was first noticed (I had looked 
in vain for possible duplicity on March 31, April 27, and June 10), 
the images in the 36-inch telescope were sharp in yellow light, being 
nearly monochromatic, and were apparently of small real diameter. 
This diameter could hardly have exceeded o%05. Further, Am was 
roughly the same for the various monochromatic emissions, indicat- 
ing that the two bodies seen had nearly identical spectra.’ This 
question was further examined in a second note,’ describing the 
visual prismatic spectrum of the two components. Van Biesbroeck’® 
has noticed the same features with the 40-inch telescope. 

Recently Aitken™ found the images “‘larger and less sharply de- 
fined when the focus was set for green rays than when it was adjusted 
for yellow light. In the latter focus the two nuclei were small and 
fairly sharp, but I should hardly call them stellar.” The seeing was 3 
during this observation. 

In view of the importance of the matter, the writer has paid spe- 
cial attention to the appearance of the images whenever a first-class 
opportunity occurred. The following records were made immediately 
after the observations with the 40-inch telescope. 

October 19, 1936.—Seeing 4. More than 95 per cent of the visual 
light is concentrated in the nebulium lines. The size of the compo- 
nents seems to be about 073 or 0’ 4 with power 960 X, and the images 
seem nearly in contact. But with 1340 X, sharp nuclei are visible— 
the true diameter of which cannot be larger than 0” 10—which are 
surrounded by hazy disks, about o”4 in diameter. The impression 
was that this larger disk (or blur) was due to the combined effect of 
imperfect correction of the 40-inch objective for blue light and of 
imperfect seeing. No continuous spectrum could be seen, or any 
nebulosity around the nova, with low power (460 X). 

March 29, 1937.—Seeing 4. The images are separated by at least 
three times their diameter. The image in the yellow plus red light 
is 4 to 5 mag. fainter than that in the “nebular” light. (Observation 


8 Pub. A.S.P., 47, 229, 1935. 10 Loc. cit. 
9 Ibid., p. 269. Loc. cit. 
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through a small prism held between the eyepiece and the eye.) The 
appearance of the yellow images, though faint, is sharp and stellar; 
the images in the nebular light look fuzzy, but that may be due to 
the objective. 

April 10, 1937.—Seeing 4. The images in the yellow light are too 
faint to be measured. The “nebular” images are fairly sharp, sepa- 
rated by at least three times their diameter. 

These observations are not conclusive as to the exact size of the 
images, but they indicate that the disks are not larger than 072, and 
are possibly only 071, in diameter. Further observations under the 
best conditions are required, and a comparison should be made with 


other objects having a similar spectrum. 
G. P. KuIPER 
YERKES OBSERVATORY 
May 1937 
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Die Fernréhre und Entfernungsmesser. By A. KOniG. 2d ed. Berlin: 

Julius Springer, 1936. Pp. v+242. Rm. 22.50. 

The second edition of this book’ contains an interesting new chapter on 
the history of the telescope from Lippershey to Porro. There is also a new 
section on instruments for measuring the altitudes of airplanes. There 
have been numerous other improvements in the text and in the illustra- 
tions. 


A Textbook of Astrophysics and Stellar Astronomy (in Russian). Edited by 
B. P. Gerasimovic. Vol. 1: Methods of Astrophysical and Astro- 
photographic Investigation.* Leningrad: ONTI, 1934. Pp. 342. Rb. 
4.20. Vol. 2: The Physics of the Solar System and Stellar Astronomy. 
Leningrad: ONTI, 1936, Pp. 579. Rb. 7.25. 

Professor Gerasimovic, who is well known in America for his brilliant 
research work in astrophysics and who is now the director of Russia’s 
largest observatory at Pulkovo, is to be congratulated for having edited 
an extraordinarily good advanced textbook of modern astrophysics. 
Russian students will have the advantage of learning from a treatise 
which contains a connected story of all phases of astrophysics. It is a pity 
that the book is not available in English. The well-known volumes by 
Russell, Dugan, and Stewart are unsurpassed as an introduction for the 
serious student; but these authors state in their Preface that ‘the pre- 
liminary knowledge assumed on the part of the student involves only the 
elements of mathematics and physics. Use of the calculus has been com- 
pletely avoided. .... ” The Russian books, on the other hand, use the 
mathematical form of exposition, although it would seem that they 
intentionally avoid the more difficult and specialized methods. 

It cannot be too strongly emphasized that our universities are par- 
ticularly in need of a textbook of the character of these volumes edited 
by Gerasimovic. Our students who specialize in astronomy rely at present 
upon Russell, Dugan, and Stewart for their general training. After that 
they begin the study of such highly specialized volumes as Eddington’s 


tAp. J., 60, 78, 1924. 


106 


| 


REVIEWS 107 


Internal Constitutions of the Stars, Rosseland’s Theoretical Astrophysics, or 
Milne’s and Pannekoek’s articles in the Handbuch der Astrophysik. The 
step from the elementary general course to the specialized advanced study 
is difficult and frequently causes even good students to lose confidence in 
their ability to master the subject. The excellent book by the two 
Strémgrens, father and son, is a useful supplement, and it includes in its 
scope not only astrophysics but also positional astronomy and celestial 
mechanics; but it is at present not available in English. 

The question arises why we have no suitable textbook which would 
serve as a link between Russell, Dugan, and Stewart and the specialized 
monographs. Is it because commercial publishers cannot make a suffi- 
cient profit from an advanced textbook in astronomy and because the 
various university presses are being operated on a profit basis? If so, we 
can perhaps make a serious effort to induce our universities, or their 
presses, to disregard the financial aspects and to publish for the few 
advanced students in astronomy the kinds of textbook they now lack. 

The volumes under review consist of a series of chapters written by 
the leading Russian astronomers in each field. 

Volume 1 begins with an Introduction by Ambarzumian on atomic 
physics. Balanovsky gives a short but complete account of the theory of 
the refractor and the reflector. Kostinsky discusses the chemical founda- 
tions of photography and contributes a long chapter on astrometric 
methods applied to photography. All necessary details are given for the 
reduction of star positions. Schlesinger’s method of dependences is also 
described. Shajn contributes a splendid chapter on the methods of stellar 
spectroscopy. A student who has carefully read this chapter should 
be able to design a good spectrograph, test it, and use it for the deter- 
mination of radial velocities. Balanovsky writes on photometry, while 
Kosirev discusses the methods and some of the results of spectrophotome- 
try. The volume is concluded by a chapter on solar work written by 
Perepelkin. 

Volume 2 contains a long chapter on the results of solar investigations 
by Perepelkin. The second chapter on the planetary system by Fessen- 
koff is especially valuable since it contains the fundamentals of theoretical 
photometry. Gerasimovic and Shajn collaborate in a chapter on “The 
Introduction to Stellar Astronomy.” This is followed by Shajn on 
“Double Stars” and by Ambarzumian on “Stellar Atmospheres.” Both 
chapters are authoritative as well as comprehensive. The remainder of 
the book is by Gerasimovi¢. It is devoted to variable stars, diffuse matter 
in interstellar space, clusters, the internal constitution of the stars, the 
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distribution of the stars in space, stellar motions, and extra-galactic 
nebulae. 

It is regrettable that such a fine book is printed on the very poor grade 
of paper which is so familiar to those who make use of the Russian 
scientific literature. The halftone illustrations are quite unsatisfactory 
and in some cases fail to show the characteristic features of the objects 


which they are supposed to illustrate. 
O. STRUVE 


Arc Spectrum of Iron from 8388 to 2242 A. By A. GATTERER and ~ 


J. Junxes. Castel Gandolfo (Italy): Specola Vaticana, 1935. Pp. 
10+ 21 illustrations. 


Spark Spectrum of Iron from \ 4650 to 2242 A. By A. GATTERER and 
J. Junxes. Castel Gandolfo (Italy): Specola Vaticana, 1935. Pp. 9+13 
illustrations. 


The first part of this publication represents a set of twenty-one 
photographic reproductions of the arc spectrum of iron. Approximately 
four thousand lines are shown and their wave-lengths, to the nearest 
o.o1 A, are indicated on the margins. The photographs were obtained 
with a prism spectrograph. The dispersion on the reproductions varies 
from 0.17 A/mm in the ultraviolet to 4.39 A/mm in the infrared. 

The second part contains thirteen photographic reproductions of the 
spectrum of the iron spark from \ 4242 to A 4650. Approximately two 
thousand lines are shown. 

The plates will be found very useful for the rapid identification of the 
comparison lines in a slit spectrogram. It is especially convenient to have 
the exact wave-lengths in international angstroms without having to 
consult the tables of Kayser and Konen. 


O. S. 
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